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Abstract 

Few  options  currently  exist  to  provide  propulsion  for  extremely  small  satellites 
due  to  design  constraints  on  power,  volume,  and  weight.  However,  future  operation  will 
require  a  capability  to  conduct  orbital  maneuvers,  momentum  dumping,  and  precision 
pointing  for  these  low  cost  satellites.  The  research  presented  here  represents  the  first 
effort  to  operate  and  quantify  the  performance  of  a  new  micro  plasma  thruster  design 
which  provides  a  novel  solution  to  these  disparate  competing  constraints.  The  thruster  in 
this  study  represents  a  deviation  from  traditional  circular  Hall  thruster  design  practices, 
because  it  eschews  a  central  magnetic  circuit,  which  results  in  nearly  parallel  electric  and 
magnetic  field  lines  within  the  thruster  discharge  channel.  This  design  decision  reduces 
thruster  complexity  and  thermal  susceptibility,  but  it  also  reduces  the  ionization 
efficiency.  The  cornerstone  of  this  study  involved  the  direct  measurement  of  thrust  in 
order  to  quantify  the  efficiency  and  specific  impulse  of  this  innovative  thruster.  The 
investigation  also  included  characterization  of  the  thruster  exhaust  plume,  voltage-current 
characteristics,  and  operating  limits.  Results  are  enumerated  and  suggestions  for 
improvement  provided. 
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PERFORMANCE  CHARACTERIZATION  OF  ANOVEE  PLASMA  THRUSTER  TO 


PROVIDE  A  REVOLUTIONARY  OPERATIONALLY  RESPONSIVE  SPACE 
CAPABILITY  WITH  MICRO-  AND  NANO-SATELLITES 

I.  Introduction 

The  extraordinary  cost  associated  with  placing  space  vehicles  into  orbit  (-'12,000 
dollars  per  kilogram)  has  served  to  justify  significant  research  and  development  in  order 
to  identify  the  highest  performing  and  most  efficient  sources  of  thrust  for  these  vehicles 
[1].  Over  the  last  50  years,  such  efforts  have  revealed  electric  propulsion  holds  the 
greatest  promise  to  fill  most  on-orbit  needs  and  has  become  the  only  viable  source  for 
interstellar  propulsion,  as  illustrated  by  Table  1.  Nonetheless,  no  single  propulsion 
technology  has  proven  equally  adept  at  meeting  the  needs  of  every  conceivable  mission. 
Instead,  space  vehicle  designers  must  weigh  factors  such  as  specific  impulse,  efficiency, 
and  system  weight  in  order  to  identify  the  best  thruster  for  a  specific  application. 

Table  1.  Typical  propulsion  performance  characteristics,  based  on  Humble  [2]. 


Technology _ Specific  Impulse  (s) _ Thrust  (N) 


Cold  Gas 

0.1  -  50 

N2 

60 

H2 

250 

Chemical 

0.01  -  12,000,000 

Liquid 

Monopropellant 

140  -  235 

Bipropellant 

320  -  460 

Solid 

260  -  300 

Hybrid 

290  -  350 

Electric 

0.0001  -  20 

Electrothermal 

500  -  1,000 

Electromagnetic 

1,000  -  7,000 

Electrostatic 

2,000  -  10,000 

1 


Problem  Statement 


Busek  Co.,  Inc.,  under  contract  with  the  Air  Foree  Institute  of  Teehnology,  has 
developed  a  unique  circular  magnet  thruster  design  whieh  underwent  preliminary 
charaeterization  as  part  of  this  study.  The  investigation  ineluded  ereation  of  voltage- 
eurrent  (V-I)  curves,  measurement  of  beam  current  using  a  Faraday  probe,  and  direct 
measurement  of  thrust  using  a  torsion  balanee.  Results  from  these  tests  suggest  this 
thruster,  named  the  Busek  BHT-20,  behaves  like  a  circular  micro-Hall  thruster  in  many 
ways,  but  has  some  peeuliarities  whieh  will  require  further  investigation  in  order  to  fully 
understand  both  their  cause  and  impact  on  overall  thruster  performance.  The  current 
iteration  of  this  design  includes  an  ionization  and  aceeleration  chamber  comprised  of  an 
opened  ended  eircular  section  with  the  walls  of  the  eylinder  being  made  of  boron  nitride 
(BN).  Surrounding  the  BN  insulator  is  a  cylindrical  0.135  Tesla  samarium  cobalt 
(SmCo)  permanent  magnet  from  Eleetron  Energy  Corp.  in  Eandisville,  PA. 

Encapsulating  the  assembly  is  a  simple  aluminum  housing,  as  illustrated  by  Eigure  1.  At 
the  entrance  to  the  ionization  ehamber  is  a  traditional  stainless  steel  anode  with  a  four- 
hole  frit  for  xenon  gas  distribution.  The  center  of  the  anode  is  covered  with  a  3.175  mm 
BN  post.  During  initial  testing,  this  thruster  was  paired  with  a  3. 175  mm  Busek  BHC- 
1500  cathode.  The  geometry  and  manufacture  of  this  device  is  straightforward  when 
eompared  to  traditional  Hall  effect  thrusters,  as  illustrated  by  Eigure  2;  however,  this  new 
coneept  seeks  to  address  several  ehallenges  relating  to  the  use  of  previous  designs  for 
very  small  satellites  including  power  eonsumption,  lifetime,  and  thruster  complexity. 
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a)  Front  view  of  BHT-20.  b)  Back  view  of  BFIT-20. 

Figure  1.  Views  of  the  BHT-20  thruster  illustrating  size  and  construction. 


Cylindrical  SmCo  Magnet 
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Boron  Nitride  Center  Post 


Boron  Nitride  Insulator 


Figure  2.  Cross-sectional  diagram  of  BHT-20  thruster  illustrating  key 
components,  taken  from  Paintal  [3]. 
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The  layout  of  this  thruster  and  the  use  of  a  permanent  magnet  are  expeeted  to 
reduce  internal  degradation,  as  compared  to  traditional  designs,  and  enable  acceptable 
efficiencies  for  the  size  class  in  which  it  operates.  The  primary  failure  mechanism  for  a 
typical  Hall  effect  thrusters  is  sputtering  erosion  caused  by  plume  impingement  on  the 
walls  and  mouth  of  thruster  that  slowly  erode  the  insulator  protecting  the  magnet 
assembly  [4].  The  strength  and  orientation  of  the  magnetic  field  created  by  the  SmCo 
magnet  is  sufficient  to  induce  a  magnetic  pressure  that  contains  the  preponderance  of  the 
plasma  within  a  narrow  channel  in  the  center  of  the  thruster,  as  illustrated  in  Figure  3. 
Initial  tests,  representing  more  than  70  hours  of  operation,  do  not  show  any  visible  signs 
of  erosion,  which  serve  to  justify  further  investigation  into  this  point.  The  use  of  a 
permanent  magnet  is  also  ideal  for  a  micro-Hall  thruster  -  reducing  the  complexity  of  the 
thruster  and  improving  the  overall  efficiency  of  the  propulsion  system  by  eliminating  the 
current  necessary  to  operate  traditional  electromagnets.  Nevertheless,  this  design  does 
pose  some  unique  challenges. 


Figure  3.  Thruster  discharge  channel  during  low  power  operation. 
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Inherent  in  the  use  of  a  permanent  magnet  is  an  aeknowledgement  of  its  operating 
limits  in  order  to  aehieve  a  desirable  lifetime  for  a  low  thrust  spaee  thruster.  Over  the  last 
decade,  work  was  performed  to  vastly  improve  the  thermal  stability  and  radiation 
resistance  of  rare  earth  permanent  magnets  for  applications  such  as  ion  thrusters  and  high 
efficiency  Stirling  Radioisotope  Generators  (SRG)  [5].  A  relatively  new  class  of  ultra- 
high  temperature  (UHT)  SmCo  magnets  has  demonstrated  minimal  (less  than  0.5%) 
magnetic  irreversible  losses  at  temperatures  of  300°C  over  a  period  of  three  years.  At 
higher  temperatures  (greater  than  550°C)  even  the  most  thermally  stable  SmCo  magnets 
will  provide  reduced  flux  densities.  An  accurate  understanding  of  the  thermal 
distribution  of  this  new  plasma  thruster  during  operation  is  therefore  a  necessary 
component  of  future  investigations.  Throughout  tests  conducted  for  this  project  a  record 
of  magnetic  field  strength  was  maintained  in  order  to  determine  if  any  degradation  in 
magnet  performance  took  place  while  testing.  During  this  inaugural  study,  no  loss  of 
magnetic  field  strength  could  be  identified  using  an  F.  W.  Bell  5180  Gaussmeter.  The 
baseline  plot  of  magnetic  field  strength  along  the  centerline  of  the  plasma  discharge 
channel  is  provided  in  Figure  4.  This  plot  depicts  typical  behavior  for  a  nearly  uniformly 
charged  hollow  magnetic  cylinder.  The  concentration  of  B-field  lines  reaches  its 
maximum  at  the  center  of  the  channel  and  this  concentration  drops  off  quadratically  in 
either  direction  from  this  point. 
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Figure  4.  Magnetic  field  strength  along  thruster  discharge  channel,  taken  from 

Paintal  [3]. 

The  unusually  small  size  and  mass  flow  rate  of  this  thruster  complicates  other 
aspects  as  well;  initial  tests  have  suggested  that  the  operation  of  this  thruster  is  intimately 
related  to  placement  of  the  external  cathode.  This  placement  is  partially  due  to  the  nature 
of  this  new  thruster,  but  geometric  constraints  of  the  current  configuration  clearly  have  an 
impact.  The  outer  diameter  of  the  current  hollow  cathode  with  enclosed  geometry  keeper 
is  significant,  as  illustrated  in  Figure  5.  The  size  of  this  cathode  not  only  limits  its  viable 
orientation,  but  it  also  would  prove  troublesome  if  trying  to  integrate  the  thruster  onto  a 
small  space  vehicle,  such  as  a  CubeSat.  The  employment  of  either  a  micro-hollow 
cathode  or  a  lanthanum  hexaboride  single  crystal  emitter  relying  on  ohmic  heating  to 
generate  electrons  will  be  explored  in  future  research.  Advantages  to  the  latter  being  a 


6 


significant  portion  of  the  xenon  propellant  will  not  be  used  by  the  eathode  itself 
(significant  loss  mechanism),  as  well  as  avoiding  any  asymmetrie  thrust  created  by  the 
eathode  itself.  A  BHT-20  speeifie  eathode  would  represent  a  signifieant  improvement, 
but  it  was  deeided  that  performance  eharaeterization  of  the  existing  system  was  more 
pressing.  The  value  of  this  thruster  eoneept  ean  only  be  aeeurately  eommunicated  by 
performanee  figures  of  merit:  thrust,  specific  impulse,  and  various  measures  of 
effieieney.  The  measurement  of  these  attributes  is  based  on  traditional  Hall  thruster 
relationships  and  theory. 


Figure  5.  BHT-20  during  operation,  illustrating  the  large  relative  eross-seetion  of  the 

eathode  to  the  plasma  plume  itself. 

Research  Objectives 

The  objeetive  of  this  researeh  effort  is  to  quantify  the  performanee  of  this  novel 
deviation  from  traditional  micro-Hall  thruster  designs,  in-terms  of  thrust,  speeifie  impulse 
(Isp),  and  effieieney.  This  researeh  will  use  xenon  for  the  propellant  in  the  BHT-20.  The 
thrust  produeed  will  be  measured  using  a  torsion  balanee  system,  whieh  ean  resolve 
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micro-Newton  forces.  The  thrust  measurement  in  conjunction  with  measurements  of 
propellant  flow,  cathode  potential,  and  anode  potential  will  be  used  to  quantify  Isp  and 
efficiency.  A  Faraday  probe  will  be  used  to  measure  the  current  density  within  the 
plasma  plume  at  various  radii  in  order  to  determine  the  ion  distribution  and  beam 
efficiency.  Erosion  of  the  boron  nitride  insulating  material  and  magnetic  field  strength 
will  be  monitored  in  order  to  determine  operational  limits  for  the  thruster  and  to 
determine  viability  for  prolonged  usage  in  the  space  environment.  This  information  will 
factor  into  a  judgment  regarding  the  viability  of  this  technology  to  fulfill  various 
applications  within  the  space  propulsion  regime  to  include  precision  attitude  control  for 
ESPA  class  and  larger  satellites  and  primary  propulsion  for  CubeSats. 


II,  Theory  and  Previous  Research 


This  chapter  provides  baekground  information  on  the  operation  of  Hall  effeet 
thrusters  and  methods  to  eharaeterize  their  performanee.  A  diseussion  of  the  physies 
governing  Hall  effeet  thrusters  is  ineluded.  Typical  applications  for  this  thruster  design 
and  inherent  limitations  are  also  explored.  The  diagnostie  teehniques  that  are  used  in  this 
study  are  introdueed  along  with  experimental  considerations.  Finally,  the  relationships 
deseribing  Hall  thruster  effieiency  are  introdueed. 

Hall  Thrusters 

Hall  thrusters  are  ineluded  under  the  broad  eategory  of  eleetrie  propulsion 
systems,  beeause  the  meehanism  by  whieh  these  devices  generate  foree  involves 
imparting  eleetrie  energy  to  ionized  partieles,  whieh  in-turn  are  aeeelerated  to  high 
veloeity  in  order  to  aehieve  an  effieient  momentum  transfer.  Electrostatie  and 
eleetromagnetie  thrusters  are  the  two  most  eommon  forms  for  an  eleetrie  propulsion 
device.  Both  approaches  generate  thrust  using  the  Lorentz  force  [2]. 

F  =  q{E  +  vxB)  (1) 

The  Lorentz  foree  is  equal  to  the  product  of  the  eharge  of  the  moving  particles  (q) 

and  the  sum  of  the  applied  eleetrie  field  (e)  and  the  eross  produet  of  the  partieles 
veloeity  (fi)  and  the  magnetie  field  strength  (b).  Hall  effeet  thrusters  and  ion  thrusters 
are  termed  eleetrostatie  thrusters  beeause  they  use  the  qE  eomponent  of  the  Lorentz  foree 
to  generate  thrust.  Hall  effeet  thrusters  transfer  momentum  by  using  the  axial  eleetrie 
field  from  the  anode  to  aecelerate  the  eharged  ions.  The  physies  model  to  explain  this 
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process  in  more  detail  is  involved,  but  the  physical  form  it  necessitates  is  not,  as 
illustrated  by  Figure  6. 


C  athode  Neutralizer 


Discharge  chamber 


Anode 
propellant 
gas  inlet 


Anode 


©  Neutrals 
O  Ions 
©  Electrons 


Electron 

current 

(ionization) 


Electron 

current 

(neutralization) 


g  t 

O 

O, 


Ion 

current 


Magnetic 

Circuit 


O  (high-energv) 

O 

Azimuthal  electron  Hall  current 
(out  of  the  paper) 


O 

O, 


O 

o 


Figure  6.  Hall  effect  thruster  cross  sectional  view,  taken  from  Warner  [6]. 
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Hall  effect  thrusters  are  commonly  composed  of  an  annular  channel  with  walls 
made  from  an  insulator,  such  as  BN,  and  an  anode  and  propellant  dispersion  medium  at 
the  inlet.  A  cathode  is  commonly  mounted  externally  near  the  exit  of  this  channel  in 
order  to  produce  electrons  to  neutralize  the  positively  charged  ions  in  the  plasma  plume. 
This  is  necessary  to  prevent  a  charge  buildup  on  the  spacecraft  and  the  resulting  back 
streaming  of  ions  that  would  net  no  thrust  [7].  The  cathode  can  also  be  mounted  down 
the  centerline  of  the  thruster,  but  this  is  less  typical.  Typically,  within  the  inner  BN 
cylinder  and  outside  the  outer  BN  cylinder  are  permanent  or  electromagnets  that  produce 
a  radial  magnetic  field  perpendicular  to  the  BN  walls.  External  to  the  thruster  is  a  power 
supply  that  links  the  anode  and  cathode  in  order  to  generate  the  axial  electric  field. 

Hall  effect  thruster  operation  starts  from  the  fraction  of  the  electron  current 
generated  by  the  cathode  that  is  used  for  ionization.  These  electrons  flow  into  the 
ionization  channel  due  to  their  attraction  to  the  anode  (150+  V).  The  electrons  cannot 
directly  migrate  to  the  anode  due  to  Hall  effect  that  is  generated  from  the  interaction  of 
the  axial  electric  field  and  radial  magnetic  field  near  the  exit  of  the  channel.  This 
interaction  produces  a  force  on  the  electrons  in  the  azimuthal  direction  that  grows  as  they 
approach  the  channel  walls  due  to  the  increased  magnetic  field  strength  [8].  The  balance 
of  these  forces  causes  the  electrons  to  circle  the  channel  in  a  closed  drift.  Further,  the 

azimuthal  velocity  (i^^)  of  these  electrons  generates  a  force  (q(v±  x  B)  j  that  causes  the 

electrons  to  gyrate  in  the  direction  of  their  azimuthal  velocity.  The  radius  of  this  gyration 
is  known  as  the  Larmour  radius  (r^)  and  is  found  using  (2)  where  m  is  the  mass  of  the 
electron  and  q  is  the  charge  of  the  electron. 
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mPx 


(2) 


The  depth  of  the  ionization  ehannel  is  sized  sueh  that  the  Lamour  radius  for  the 
electrons  is  substantially  smaller  than  the  depth  of  the  channel  and  the  Lamour  radius  of 
the  ions  is  substantially  larger  that  the  depth  of  the  channel  [8].  Sizing  with  respect  to  the 
electron  is  done  to  produce  a  high  statistical  likelihood  of  an  ionizing  collision  between 
the  neutrals  and  the  electrons.  Sizing  with  respect  to  the  ions  is  done  to  ensure  that  the 
magnetic  field  does  not  drastically  perturb  the  direction  of  ions  acceleration  out  of  the 
channel.  The  electrons  will  eventually  work  their  way  to  the  anode  after  collisions  with 
particles  and  the  channel  walls  where  they  will  be  recollected  and  serve  to  maintain  the 
anode  current. 


The  propellant  is  fed  into  the  ionization  channel  near  the  anode  after  passing 
through  a  dispersion  media  [8].  A  large  percentage  of  the  neutral  propellant  atoms, 
xenon  being  the  most  common  choice  due  to  its  high  atomic  weight  and  low  ionization 
potential,  then  interact  with  the  electron  cloud  described  above.  The  collisions  with  high 
energy  electrons  causes  ionization  of  the  propellant,  and  these  ions  are  accelerated  by  the 
electric  field  to  high  axial  velocities  (>10  km/s).  The  ion  plume  outside  the  thruster  is 
then  neutralized  by  electrons  being  generated  by  the  cathode  in  order  to  preserve  a  quasi¬ 
neutral  plasma.  The  theory  of  operation  for  a  traditional  Hall  thruster  does  not 
completely  capture  the  physics  underlying  the  BHT-20;  nonetheless,  all  of  the 
performance  measurements  developed  from  the  theory  outlined  above  are  still  valid  for 
this  thruster.  The  ability  to  make  this  assumption  is  important,  because  it  allows  this 
investigation  to  build  on  nearly  50  years  of  Hall  effect  research. 
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Hall  effect  thruster  development  began  in  1960’s  in  both  the  United  States  (US) 
and  the  former  Soviet  Union  [9].  The  United  States  then  abandoned  research  involving 
this  thruster  for  more  than  two  decades  after  tests  in  the  late  1960’s  suggested 
significantly  higher  values  for  hp  could  be  generated  by  continued  development  of 
gridded  ion  engines.  Alternatively,  researchers  in  the  Soviet  Union  elected  to  continue  to 
research  Hall  thruster  physics  and  subsequently  made  improvement  in  both  thrust  and 
efficiency.  The  Union  of  Soviet  Socialist  Republics  (USSR)  continued  this  development 
to  the  point  where  one  manufacturer,  the  Experimental  Design  Bureau/Fakel,  has 
designed  and  manufactured  more  than  110  flight  Hall  thruster  [10].  A  majority  of  these 
thrusters  were  used  for  north-south  station-keeping  (NSSK)  missions  with  one  variant  of 
the  SPT-100  demonstrating  more  than  7,000  hours  of  operation  in  orbit.  The  United 
States  eventually  become  aware  of  these  successes  and  interest  developed  in  reevaluating 
the  Hall  effect  thruster. 

In  the  1990’s,  following  the  end  of  the  Cold  War,  institutions  within  the  United 
States  were  able  to  acquire  Hall  thrusters  from  Russia  for  evaluation  [9].  The  evaluation 
was  intended  to  validate  performance  claims,  because  the  rate  of  technology  transfer 
through  results  of  previous  research  studies  performed  in  the  former  Soviet  Union  was 
slowed  by  the  lack  of  translated  manuscripts  [11].  Thus,  the  West  needed  to  commission 
many  studies  in  order  to  understand  the  physics  and  design  of  Hall  effect  thrusters  and 
this  blossomed  into  independently  designed  US  and  European  thrusters  with  favorable 
performance.  Recent  work  conducted  by  the  Air  Force  Research  Eaboratory  (AFRL)  has 
demonstrated  that  the  latest  generation  of  Hall  effect  thrusters  are  capable  of  performing 
orbit  transfer  missions  faster  and  deliver  a  greater  payload  than  ion  engines  as  a  result  of 
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their  greater  specific  power  [9].  Subsequently,  the  use  of  such  thrusters  to  provide  the 
transfer  function  for  communication  satellite  (COMSAT)  missions  has  been 
demonstrated  [12],  The  ability  of  the  Hall  thruster  to  provide  orbit  maintenance  and 
momentum  dumping  after  the  transfer  is  an  additional  benefit.  These  successes  have 
revitalized  interest  and  investment  into  the  further  advancement  of  this  technology. 

The  realization  that  Hall  thrusters  possessed  further  untapped  potential  has 
emboldened  the  National  Aeronautics  and  Space  Administration  (NASA)  to  invest  in  the 
development  of  the  High  Voltage  Hall  Accelerator  (HIVHAC)  with  improved  fip, 
throttle-ability,  and  endurance  in  order  to  provide  an  alternative  for  cost-capped  deep 
space  science  missions  [13].  The  objective  of  this  endeavor  is  to  build  upon  the 
successes  of  two  previous  Earth  orbiting  missions  that  utilized  Hall  effect  thrusters:  the 
United  States  National  Reconnaissance  Office’s  (NRO)  Science  Technology  Experiments 
(STEX)  spacecraft  that  was  launched  in  1998  and  the  European  Space  Agency’s  (ESA) 
Small  Missions  for  Advanced  Research  in  Technology  (SMART-1)  Eunar  probe  that  was 
launched  in  2003  [14].  These  two  spacecraft  are  illustrated  in  Eigure  7.  The  ultimate 
objective  is  to  make  Hall  thrusters  a  viable  candidate  for  New  Millennium  missions 
where  the  objective  is  utilize  smaller,  lighter,  and  less  expensive  systems  to  accomplish 
scientific  missions  [15].  At  the  other  end  of  the  spectrum,  interest  is  also  growing  to 
develop  small  low  power  Hall  effect  thrusters. 

Tsuyohito  Ito  and  coworkers,  from  Princeton,  suggests  there  is  increasing 
awareness  of  the  potential  to  develop  a  relatively  low  powered  Hall  effect  thruster  (10- 
lOOW)  in  order  to  provide  a  viable  propulsion  alternative  for  small  power-limited 
satellites  [16].  Research  at  Stanford  also  substantiates  this  point  [17].  The  Russian  SPT- 
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25  Hall  effect  thruster  is  one  attempt  to  fill  such  a  void.  The  SPT-25  possesses  the 
following  performance  characteristics:  power  range  of  50-150W,  total  efficiency  of  14% 
at  60W,  and  a  specific  impulse  of  550s  at  60W  [18].  In  the  United  States,  universities 
such  as  Stanford,  MIT,  and  Princeton  have  experimented  with  variations  on  a  traditional 
annular  Hall  effect  thruster  design  in  order  to  demonstrate  thrusters  operated  at  similar 
power  levels  [19].  Stanford  used  the  geometric  principles  governing  annular  Hall  effect 
thruster  operation  and  scaled  them  appropriately,  so  their  thruster  could  operate  as  low  as 
20  W  [20].  Alternatively,  Princeton  explored  the  advantages  of  using  a  circular  Hall 
thruster  (CHT)  design  for  operation  below  100  W  [21].  It  is  plausible  that  some  of  their 
future  designs  may  be  able  to  operate  at  even  lower  power  levels  with  some  sacrifice  to 
Isp  and  efficiency.  It  would  be  feasible  to  integrate  such  a  thruster,  operating  in  the  2.5  to 
5  W  range,  onto  a  standard  3U  CubeSat. 


(a)  STEX  spacecraft,  image:  NRO  [22]  (b)  SMART-1  spacecraft,  image:  ESA  [23] 

Figure  7.  NRO  STEX  spacecraft  during  assembly  and  artist  rendition  of  the  ESA 
SMART- 1  spacecraft  on  route  to  the  moon. 

Work  in  electric  propulsion  is  an  important  enabling  technology  to  expand  the 
usefulness  of  CubeSats  by  means  of  increasing  their  achievable  AV,  providing  additional 
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options  for  their  attitude  control,  and  supplying  an  efficient  momentum  dumping 
capability  [19].  Micro  Hall  effect  thrusters  are  an  important  part  of  the  effort  to  realize 
these  goals,  because  unlike  some  well-developed  technologies,  such  as  solid  fueled 
rocket  motors,  an  electric  propulsion  system  could  simultaneously  contribute  to  all  three 
objectives.  The  development  of  such  as  system  would  produce  a  paradigm  shift  in  the 
constraints  surrounding  what  is  currently  feasible  with  these  flexible  low  cost 
experimental  satellites.  However,  several  challenges  must  be  surmounted  before  the 
benefits  afforded  by  electric  propulsion  (EP)  are  realized  for  these  small  satellites. 

Among  the  engineering  trades  that  must  be  made  when  optimizing  electric 
thrusters  is  the  amount  of  thrust  that  can  be  generated  versus  the  electric  power  that  can 
be  allocated  to  the  propulsion  system  [24].  The  limited  availability  of  electric  power 
ultimately  limits  the  achievable  thrust  levels  for  these  devices.  This  is  especially  true  on 
small  power  limited  satellites,  such  as  CubeSats.  The  efficiency  of  a  Hall  effect  thruster 
tends  to  degrade  exponentially  with  decreased  size.  Efficiencies  for  thrusters  such  as  the 
Busik  BHT-200  at  200  W  is  around  42%,  whereas  the  efficiency  for  SPT-25  identified 
above  was  14%.  However,  even  lower  electrical  efficiencies  may  be  acceptable  for  small 
satellite  applications  where  the  non-electric  propulsion  alternatives  include  cold  gas 
thruster  and  butane  micro-thrusters  with  fip  values  below  70  s  [19]. 

A  consequence  of  a  Hall  effect  thruster’s  high  fip  and  low  thrust  relative  to 
traditional  chemical  thrusters  is  the  Hall  thruster  must  operate  for  significantly  longer 
durations  in  order  to  produce  the  same  change  in  velocity  (dP).  This  implies  Hall 
thrusters  require  a  longer  timeframe  to  achieve  a  given  AV  and  must  be  designed  for 
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longer  lifetimes  than  traditional  thrusters.  In  order  to  better  understand  the  principle 
failure  mode  affecting  these  thrusters  researchers  at  the  University  of  Michigan  have  been 
working  to  build  a  hydrodynamic  plasma  model  for  simulating  sputtering  damage  [25]. 
This  analysis,  as  well  as,  research  performed  at  the  NASA  Glenn  Research  Center  (GRC) 
has  demonstrated  that  the  erosion  profiles,  in  terms  of  thruster  characteristics,  are 
primarily  a  function  of  discharge  voltage,  magnetic  field  density,  and  magnetic  field 
strength  [26].  Erosion  rate,  in  terms  of  propellant  characteristics,  is  primarily  a  function 
of  ion  mass,  atomic  number,  and  the  ionization  ratio  [27].  The  propellant  for  the  Busik 
thruster  being  investigated  in  this  study  is  fixed;  however,  it  will  be  important  to 
characterize  the  thruster  parameters  and  measure  erosion  as  a  function  of  time  for  the 
thruster  in  order  to  understand  how  the  nontraditional  magnetic  field  lines  generated  by 
the  permanent  magnet  influences  thruster  lifetime.  Specifically,  it  will  be  interesting  to 
determine  if  the  increased  ion  current  density  will  significantly  shorten  this  thruster’s  life 
as  it  does  with  other  small  thruster  designs  [28]. 

A  consideration  that  must  be  taken  into  account  when  designing  a  spacecraft  to 
include  Hall  effect  thruster  propulsion  is  thoughtful  integration  to  prevent  high  energy 
divergent  ions  from  damaging  the  spacecraft’s  components,  such  as  solar  panels  [29]. 
Divergent  particles  near  the  edge  of  the  plume  can  erode  surface  finishes  such  as  those  on 
solar  panels,  reducing  power  absorption  if  these  parts  are  within  the  path  of  the  particles 
[30].  Surfaces  that  are  treated  using  electrostatic  discharge  (ESD)  techniques  are  also 
susceptible.  The  size  of  this  plume  also  has  implications  for  satellite  operating 
procedures.  The  distribution  of  charged  particles  can  slightly  increase  the  signal  to  noise 
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ratio  of  communication  with  ground  stations  while  the  thruster  is  operating.  In  both  oases 
the  obstaoles  are  surmountable  with  prior  planning. 

Variants  on  the  Annular  Hall  Thruster 

It  is  relatively  diffioult  to  soale  annular  Hall  thruster  designs  down  to  low  power 
due  to  the  need  to  hold  the  main  dimensionless  parameters  oonstant.  This  implies 
deoreasing  the  ohannel  size  and  inoreasing  the  magnetio  field  strength.  In-tum,  this 
signifioantly  oomplioates  the  design  of  the  magnetio  oirouit,  reduoes  the  volume  available 
to  integrate  heat  shields  and  erosion  resistant  insulation,  and  oan  lead  to  signifioant  heat 
loading  and  erosion  of  ohannel  surfaoes  [31].  The  one  notable  suooessful  attempt  to  soale 
an  annular  Hall  effeot  thruster  to  operate  at  less  than  50  W  was  done  by  Stanford 
University  [32].  This  design  inoluded  diamond  lined  ohannel  walls  in  order  to  provide 
improved  temperature  and  sputtering  erosion  resistanoe.  Previous  researoh  performed  at 
Stanford  on  an  inherently  less  effioient  linear  geometry  Hall  thruster  had  demonstrated 
that  diamond  walls  provided  a  25%  deorease  in  the  sputtering  erosion  rate  when 
oompared  with  traditional  boron  nitride  walls  [33].  The  one  negative  attribute  assooiated 
with  using  diamond  was  the  extremely  high  thermal  oonduotivity.  This  is  one  of  the 
reasons  the  Stanford  annular  design  had  to  inelude  signifioant  water  oooling  of  the 
thruster  magnets  in  order  to  avoid  the  demagnetization  that  would  have  resulted  from  the 
unoooled  operating  temperature.  These  ohallenges  led  researohers  at  the  Prinoeton 
Plasma  Physios  Laboratory  to  pursue  development  of  a  low  power  oiroular  Hall  thruster. 

A  oylindrioal  Hall  thruster  has  a  single  oentral  magnet  that  is  generally  flush  or 
nearly  flush  with  the  anode,  whioh  deoreases  the  insulator  surfaoe  area  to  disoharge 
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channel  volume  and  simplifies  eonstruction  for  low  power  applieations  [34],  Similar  to 
annular  Hall  thrusters,  CHTs  may  use  either  permanent  or  electromagnets  to  generate 
quasi-orthogonal  magnetic  field  lines  between  the  central  circuit  and  outer  discharge 
ehannel  wall.  The  axial  component  of  magnetie  field  lines  is  signifieantly  more 
pronounced  than  for  annular  designs,  but  the  general  physies  underlying  CHT  operation 
does  not  appreciably  change  [35].  The  BHT-20  ean  be  seen  as  a  further  evolution  of  the 
CHT  design  where  the  eentral  magnetic  circuit  is  removed  to  ease  packaging  as  the 
thruster  is  scaled  to  very  low  power  levels. 

The  fundamental  eleetromagnetic  differenee  between  the  BHT-20  and  a 
traditional  CHT  is  the  magnetic  field  lines  for  the  BHT-20  are  nearly  axial  in  orientation. 
The  lack  of  a  significant  radial  component  to  the  magnetic  field  lines  implies  that  no 
significant  Hall  effect  will  be  produced  by  eleetrons  following  the  primarily  axial  eleetrie 
field  lines  to  the  anode.  It  is  supposed  that  this  will  reduce  the  electron  dwell  time  in  the 
discharge  ehannel  and  possibly  reduee  the  beam  current  efficiency  as  a  result. 

Diagnostic  Techniques 

Diagnostic  instruments  identified  for  this  researeh  include  two  different  thrust 
stand  designs  and  a  Faraday  probe.  These  diagnostie  tools  are  well  established  within  the 
field  of  eleetrie  propulsion  and  have  been  used  extensively  at  AFIT  to  characterize  both 
established  and  new  thruster  designs.  The  use  of  well  vetted  instruments  was  deemed 
important  due  to  the  laek  of  previous  knowledge  regarding  the  operation  of  this  novel 
micro  thruster  design.  The  Faraday  probe  was  the  first  instrument  used,  beeause  the 
laboratory  was  already  configured  for  this  apparatus. 
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At  its  core,  the  Faraday  probe  consists  of  a  metal  collector  plate  that  measures  ion 
current  at  various  points  in  the  plasma  plume.  The  aceelerated  ions  in  the  plume  strike  a 
collector  plate  and  this  induces  a  flow  of  electrons,  which  is  the  eurrent  that  the  probe 
measures.  In  order  to  reduee  one  major  souree  of  error,  the  Faraday  probe  used  in  this 
experiment  ineludes  a  guard  that  can  be  biased  negative  in  order  repel  the  electrons  in  the 
plume,  and  the  ring  reduces  distortion  of  the  sheath  by  the  probe.  These  eleetrons  would 
artificially  reduce  the  measured  eurrent  if  they  were  to  strike  the  eollector  plate.  The 
measured  eurrent  and  the  guard  port  area  ean  be  used  to  ascertain  the  eurrent  density  at  a 
point  in  the  plasma  plume.  A  map  of  the  current  density  in  a  plane  ean  be  ereated  once  a 
series  of  measurements  are  taken  at  a  eonstant  radius  from  the  eenterline  of  the  thruster 
discharge  ehannel.  The  result  provides  insight  into  plume  symmetry,  divergence  angle, 
and  with  numerie  integration,  the  beam  eurrent.  All  of  these  performanee  measures  help 
to  reveal  how  the  thruster  funetions  and  how  well  it  funetions,  but  a  measurement  of 
thrust  is  still  neeessary  in  order  to  identify  the  fundamental  performanee  parameters  for  a 
thruster. 

The  most  important  performance  measurement  that  is  used  to  eharaeterize  a  new 
propulsion  deviee  is  thrust.  In  this  experiment,  an  attempt  was  made  to  measure  thrust 
with  two  different  instruments.  The  average  propellant  veloeity  (v)  and  Igp  can  be 
directly  derived  using  the  measured  value  for  thrust  (T)  and  the  mass  flow  rate  of  the 
propellant  (rh)  [36]. 
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The  thrust  efficiency  can  also  be  determined  using  T  once  the  total  system  power  P  is 


measured. 
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The  first  instrument  that  was  used  to  measure  the  thrust  output  of  the  BHT-20  was 
a  traditional  inverted  pendulum  thrust  stand.  This  stand  design  is  commonly  used  to 
characterize  low  thrust  (mN)  electric  propulsion  devices,  such  as  Hall  effect  thrusters. 
These  instruments  can  generally  be  tuned  to  accommodate  a  wide  variety  of  thruster 
weights  and  outputs  through  the  appropriate  selection  of  eight  flexures  and  a  primary 
spring.  The  accuracy  of  this  stand  is  highly  dependent  on  the  validity  of  an  assumption 
that  all  measurements  are  taken  at  a  constant  temperature.  The  use  of  active  water 
cooling  and  frequent  calibrations  is  the  most  common  way  to  reduce  this  potential  source 
of  uncertainty  [37]. 

The  second  instrument  that  was  used  to  measure  the  thrust  output  of  the  BHT-20 
was  a  torsional  thrust  balance  built  by  Busek,  which  is  capable  of  resolving  continuous 
thrust  levels  of  approximately  IpN  [38].  This  design  requires  the  thruster  to  be  mounted 
on  one  side  of  a  balance  arm  with  an  appropriate  counterweight  on  the  opposite  side  in 
order  to  balance  the  arm  and  minimize  any  torques  generated  by  the  weight  of  the 
thruster.  Motion  of  the  balance  arm  is  resisted  by  a  calibrated  torsional  spring,  and  the 
displacement  is  measured  using  a  laser  displacement  system  (LDS).  Within  the  operating 
range  of  the  torsional  balance,  the  LDS  is  capable  of  producing  a  nearly  linear  voltage 
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response  curve  with  respect  to  displacement.  A  schematic  of  the  torsional  thruster  stand 
and  the  associated  support  equipment  is  illustrated  in  Figure  8. 


Figure  8.  Schematic  of  torsion  thrust  balance  for  measuring  pN  thrust  levels. 

The  torsional  balance  provides  a  precise  technique  to  measure  thrust,  but  the 
accuracy  of  the  derived  thrust  values  are  highly  sensitive  to  the  accuracy  of  the 
calibration  voltage  that  is  generated  during  setup.  The  voltage  reading  will  drift  over 
time  if  the  system  does  not  operate  at  a  constant  temperature.  This  necessitates  checking 
the  system  zero  before  and  after  each  measurement  in  order  to  determine  if  any  drift 
occurred.  If  drift  was  present  then  a  linear  approximation  of  the  drift  rate  provides  a 
simplified  method  to  account  for  the  influence  on  the  reported  voltage  values  with  respect 
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to  time.  Thrust  measurements  at  eonstant  power  simplify  the  analysis  by  avoiding 
difficulties  arising  from  resolving  transient  loads  with  this  instrument. 

Experimental  Considerations 

AFIT  currently  possesses  a  3.175  mm  Busek  BHC-1500  hollow  cathode  with 
enclosed  keeper.  While  testing  at  Busek  has  demonstrated  that  this  cathode  will  work 
with  the  thruster  under  investigation,  the  packaging,  propellant  utilization  efficiency,  and 
thrust  vectoring  that  result  are  not  desirable.  One  alternative  would  be  to  acquire  a 
micro-hollow  cathode  of  similar  design,  but  with  flow  channel  diameter  of  1.5875  mm  or 
smaller.  The  cost  associated  with  this  option  is  prohibitively  expensive  and  the  timeline 
associated  with  development  of  such  a  device  is  currently  unknown.  The  other 
alternative  appears  more  feasible  at  this  point  in  time. 

The  hollow  cathode  could  be  replaced  by  a  large  lanthanum  hexaboride  single 
crystal  emitter  relying  on  ohmic  heating  to  generate  electrons.  Lanthanum  hexaboride 
(LaBe)  Cathodes  offer  the  capability  of  long  life  relative  to  alternative  electron  emitters 
and  are  more  than  two  orders  of  magnitude  less  sensitive  to  impurities  present  in  the 
propellant  than  traditional  alternatives  [39].  The  longevity  of  LaBg  has  been  confirmed 
by  tests  including  operation  for  over  300  h  while  being  exposed  to  air  more  than  100 
times  without  measureable  degradation  in  emission  [40].  Figure  9  offers  a  view  of 
multiple  LaBe  emitters.  The  emitter  pictured  in  the  lower  left-hand  comer  of  the  image  is 
suitable  for  use  with  low  current  space  propulsion  devices  due  to  the  relatively  large 
surface  area  of  the  crystal  that  can  be  ablated  to  produce  electrons. 
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Figure  9.  Depiction  of  multiple  Lanthanum  Flexaboride  crystal  emitter  designs, 
used  with  permission  of  Applied  Physics  Technologies. 

Efficiency  Determinations 

The  total  efficiency  of  a  Hall  effect  thruster  was  previously  described  using  the  jet 
power  relationship.  The  total  efficiency  can  also  be  expressed  as  the  decomposition  of 
the  component  efficiencies,  as  shown  in  (6)  from  Goebel  and  Katz  [8]. 

Vt  =  Y^VbVvVmVo  (6) 

The  correction  factor  (y)  is  used  to  describe  the  effects  of  beam  divergence  and 
multiply  charged  ion  content  and  is  given  by 

Y^ocFf  (7) 
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The  thrust  correction  factor  (a)  is  used  to  account  for  the  presence  of  doubly  ionized 
atoms.  Letting  represent  the  number  of  singly  ionized  atoms  and  represent  the 
number  of  doubly  ionized  atoms  produces  the  following  expression  for  a, 


a  = 


1+ 


/++ 
V2/  + 


1  + 


/+  + 
/  + 


(8) 


Thus,  the  ratio  of  singly  to  doubly  ionized  atoms  can  be  found  directly  using  laser 
absorption.  The  geometric  efficiency  (F(-)  represents  the  thrust  loss  from  beam 
divergence  and  is  found  using  the  half-angle  divergence  of  6  and  the  ion  current  density 
(J (r)),  as  shown  in  (9). 

C  2nJ(r)  cos  6  (r)dr 

Ft  =  ^  (9) 

Ib 


The  /(r)  distribution  can  be  found  experimentally  using  a  Faraday  probe  at  varying 
values  of  r.  The  beam  current  efficiency  (r;^)  is  found  from  the  ratio  of  the  beam  current 
to  the  discharge  current. 

Vb  =  r  (10) 

Id 


The  voltage  ratio  (j]p)  is  found  from  the  ratio  of  the  beam  voltage  to  the  discharge 
voltage. 

(11) 
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The  mass  utilization  efficiency  (r/m)  is  found  from  the  ratio  of  the  mass  flow  of  ions  out 
of  the  thruster  (rhj)  to  the  sum  of  the  mass  flow  rates  from  the  anode  (rh^)  and  cathode 

(rhc). 
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The  electrical  utilization  efficiency  (rjo)  is  found  from  the  ratio  of  the  discharge  power  of 
the  thruster  (P^)  to  the  sum  of  P^,  the  keeper  power  (P^),  and  magnet  power  {Pmag}- 
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The  total  efficiency  of  this  p,-Hall  thruster  will  be  found  with  the  jet  power 
relationship.  The  other  isolated  efficiencies  above  will  be  reported  based  on 
available  information  from  testing  in  order  to  better  characterize  this  thruster. 
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Ill,  Methodology 


This  chapter  details  the  methodology  and  experimental  setup  used  to  generate 
data  for  this  new  thruster  design.  Details  regarding  the  laboratory  equipment  and  faeility 
used  to  conduct  this  research  are  provided.  The  proeedures  used  to  setup  and  eonduet 
eaeh  experiment  are  delineated.  The  methods  used  to  reduce  the  data  from  eaeh 
experiment  and  the  methods  used  to  quantify  the  experimental  error  are  deseribed. 

Vacuum  Facility  and  Support  Equipment 

The  experimental  work  detailed  within  this  paper  took  place  within  the  AFIT 
Spaee  Propulsion  Analysis  and  System  Simulation  (SPASS)  Laboratory.  This  faeility 
eontained  the  largest  and  most  versatile  vacuum  chamber  at  AFIT.  The  experimental 
setup  eentered  around  this  2.5  meter  long  by  1.8  meter  vaeuum  ehamber  that  was  eapable 
of  aehieving  vacuum  levels  on  the  order  of  10'^  torr.  The  chamber  was  produced  and 
installed  by  PHPK  Teehnologies.  It  featured  multiple  feed  through  and  view  ports  that 
made  it  an  ideal  facility  for  supporting  the  instrumentation  used  to  eonduet  this  test,  as 
illustrated  in  Figure  10.  In  order  to  expedite  deereasing  the  tank  pressure  from 
atmospherie  to  10'^  torr  a  staged  system  was  utilized. 

The  first  phase  of  the  pump  down  sequenee  relied  on  a  meehanieal  ruffing  pump, 
a  Leybold  Trivae  D65B,  to  reduee  pressure  within  the  ehamber  from  atmospheric  to  its 
crossover  pressure  (100  millitorr)  at  whieh  point  it  would  deaetivate.  Upon  initiation  of 
the  pump  down  sequenee,  a  helium  eooled  CVI  Torr  Maser  TM250  began  eooling.  This 
eryogenie  pump  was  mounted  atop  the  ehamber,  but  separated  from  the  tank  with  a  gate 
valve  that  would  not  open  until  the  erossover  pressure  was  aehieved.  The  TM250  then 
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began  reducing  the  pressure  within  the  chamber  farther  while  four  helium  cooled  CVI 
Torrmaster  TM500  cryo-pumps  began  cooling.  Upon  reaching  their  operating 
temperature,  which  ranged  from  15  to  20  Kelvin,  the  combined  pumping  capacity  of 
these  four  pumps  was  16,000  liters  per  second  of  xenon  gas  [38].  This  translated  to  a 
typical  operating  pressure  of  approximately  10'^  torr  while  conducting  experiments  with 
the  BHT-20.  Accurate  measurement  of  the  pressure  during  operation  was  essential  to 
ensure  inadvertent  cathode  poisoning  did  not  occur  and  to  be  able  to  record  the  pressure 
when  experiments  were  being  conducted. 


Figure  10.  AFIT  SPASS  laboratory  chamber  and  support  equipment. 


Two  devices  were  used  to  monitor  pressure  and  measure  the  gas  species  present 
in  the  vacuum  chamber  during  testing.  A  Lesker  300  Series  Convection  vacuum  gauge 
was  used  to  measure  pressures  in  excess  of  0.1  millitorr,  which  was  insightful  during 
roughing  pump  operation.  An  Extorr  XT  100  Residual  Gas  Analyzer  (RGA)  was  used  to 
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5  8 

measure  pressures  from  10'  to  10'  torr.  This  unit  integrated  three  separate  gauges,  so  it 
was  capable  of  measuring  pressures  ranging  from  atmospheric  to  10'^^  torr.  Extorr’s 
VacuumPlus  software  was  used  during  testing  to  provide  both  a  continuous  pressure 
reading  and  a  plot  of  the  partial  pressure  for  the  gas  species  present  in  the  vacuum 
chamber,  as  shown  in  Figure  1 1 .  Special  care  was  taken  to  ensure  H2O  and  O2  levels 
were  below  10'^  before  cathode  conditioning  commenced  to  avoid  uimecessary  damage. 
The  automation  built  into  this  system  helped  to  decrease  the  opportunity  for  human  error, 
increased  test  reproducibility,  and  expedited  the  pumping  sequence. 
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Figure  1 1 .  Screenshot  of  the  Extorr  VacuumPlus  software  interface  with 
identification  of  significant  gaseous  species  components. 
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The  vacuum  chamber  was  regularly  pumped  up  to  atmospheric  pressure  after  the 
completion  of  a  test  in  order  to  check  the  condition  of  the  thruster’s  boron  nitride 
insulator  and  permanent  magnet.  The  procedure  was  not  initiated  until  the  cathode  had 
ceased  operation  for  at  least  one  hour  in  order  to  ensure  that  the  emitter  was  not  hot 
during  the  pump-up  sequence.  To  initiate  the  sequence  a  single  stop  button  was 
depressed  that  caused  the  five  cryo-pumps  to  cease  operation.  The  chamber  was  filled 
with  dry  nitrogen  once  the  CVI  CBST  6.0  scroll  compressors  ceased  cooling  the  pumps 
and  the  gate  valve  for  the  TM250  closed.  The  dry  nitrogen  was  supplied  via  a 
compressor  and  battery  of  nitrogen  generators.  Once  the  chamber  was  at  atmospheric 
pressure  the  battery  of  instruments  and  tools  within  the  chamber  were  easily  accessible 
for  establishing  new  experimental  setups. 

To  facilitate  the  precise  placement  of  the  cathode  relative  to  the  thruster  and  other 
instruments  within  the  chamber,  such  as  the  Faraday  probe,  the  cathode  was  mounted  on 
a  fixture  constructed  of  80/20'^’’^  aluminum  components  that  were  generally  bolted  to  the 
vertical  Z-axis  of  a  three-axis  Aerotech®  translation  system.  The  thruster  was  oriented 
pointing  down  the  length  of  the  chamber  and  it  was  fixed  in  this  orientation,  because  the 
translation  system  could  only  move  in  the  three  primary  axes  which  were  oriented  such 
that  the  X-axis  was  parallel  to  the  length  of  the  chamber.  Each  axis  of  the  translation 
stage  possessed  a  60  cm  range  of  motion  with  sub-millimeter  accuracy  on  the  movement 
of  each  axis.  The  flexibility  afforded  by  this  system  made  it  possible  to  quickly  transition 
the  cathode  from  one  position  to  another  relative  to  the  thruster  without  opening  the 
chamber,  as  depicted  in  Figure  12.  Because  the  thruster  position  was  fixed,  the  Faraday 
probe  also  had  to  be  mounted  to  a  translation  system,  as  illustrated  by  Figure  13. 
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a)  Side  view  of  cathode  translation  mount  b)  Front  view  of  cathode  translation  mount 
Figure  12.  Cathode  mounted  to  translation  stage  for  movement  relative  to  thuster. 


Figure  13.  Faraday  probe  translation  stage  depicted  in  foreground. 

The  Faraday  probe  was  mounted  on  a  three  axis  translation  and  rotation  system, 
which  was  designed  by  Colorado  State  University.  The  system  included  a  radial  stage 
(R),  an  axial  stage  (Z),  and  a  rotational  stage  (0).  The  R  and  Z  stage  both  possessed  a  1 .0 
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meter  range  of  motion  and  the  0  stage  possessed  a  180°  range  of  motion.  The  three 
stepper  motors  in  this  system  were  water  cooled  by  a  Thermo  Scientific  NESLAB  RTE7 
chiller  external  to  the  vacuum  chamber.  The  stages  and  Earaday  probe  were  all 
controlled  by  a  single  integrated  EabView®  VI  developed  by  Colorado  State  University 
for  this  purpose. 

BHT-20  Control  System 

The  BHT-20  thruster  was  operated  over  a  range  of  voltages  and  discharge 
currents  that  were  found  to  be  stable  during  the  generation  of  voltage-current  curves.  The 
thruster  was  paired  with  an  external  BHC-1500  cathode  that  was  operated  at  the 
manufacturer’s  suggested  parameters.  An  individual  xenon  propellant  feed  line  was  used 
for  both  the  thruster  and  cathode  during  operation,  and  both  devices  were  controlled 
using  a  Busek  Power  Processing  Unit  (PPU).  The  operating  procedures  established  by 
the  author  for  the  BHT-20  are  included  in  Appendix  A. 

The  SPASS  laboratory  was  plumbed  to  provide  either  xenon  or  krypton  for  use  in 
electric  thrusters.  The  xenon  gas  used  in  this  experiment  was  stored  in  a  pressurized  tank 
with  a  regulator  that  was  used  to  reduce  the  line  pressure  from  approximately  600  psi  to 
20  psi.  The  xenon  propellant,  flowing  at  20  psi,  was  then  routed  to  four  MKS  model 
180A  precision  mass  flow  controllers,  as  shown  in  Eigure  14.  Two  of  these  mass  flow 
controllers  had  a  full  range  of  50  seem  and  two  had  a  range  of  10  seem.  All  of  the 
controllers  have  an  accuracy  of  0.01%  of  their  full  flow  value  [41].  Both  the  cathode  and 
thruster  used  a  10  seem  controller.  Xenon  flow  was  routed  from  these  controllers  to  the 
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thruster  and  cathode  with  a  combination  of  stainless  steel  and  PTFE  lines.  Flow  rates 
were  controlled  and  monitored  using  an  MKS  Type  247  four-channel  readout. 


Figure  14.  Xenon  and  krypton  bottle  and  battery  of  four  mass  flow  controllers. 


Power  for  the  thruster  and  cathode  were  provided  by  a  Busek  BPFl-600  Power 
Processing  Flnit  (PPFl),  which  reported  discharge  voltage  with  an  accuracy  of  0.1%  [42]. 
The  PPU  received  its  power  from  a  Sorensen  DCS55-55  power  supply,  which  was 
capable  of  0-55  V  and  0-55  A.  Control  of  the  PPU  was  achieved  using  Busek’s  BPU-600 
Host  Simulator,  which  was  a  FabView®  software  interface  and  is  shown  in  Figure  15. 
This  program  provided  both  the  capability  to  control  the  functions  associated  with  the 
thruster  and  cathode,  as  well  as,  monitor  their  performance. 
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Figure  15.  BPU-600  Flost  Simulator  software  interface. 


Voltage-Current  Measurement 

In  order  to  generate  voltage-current  (V-I)  curves,  the  voltage  and  current  for  the 
thruster  and  cathode  had  to  be  measured  and  recorded.  The  anode  voltage  was 
commanded  and  monitored  using  Busek’s  BPU  Host  Simulator  Software.  Shunt  resistors 
connected  in  series  were  used  in  conjunction  with  Ohm’s  Law  to  measure  current  for  the 
anode  and  cathode  keeper.  A  National  Instruments  SCXl-1321  4-channel  readout  was 
used  to  measure  the  voltage  drops  across  the  resistors,  and  these  voltages  were  used  to 
generate  current  values  within  the  LabView®  program  for  controlling  the  PPU.  A 
special  80/20™  mounting  structure  for  the  thruster  and  cathode  was  produced  that 
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enabled  the  thruster  to  face  one  portal  while  it  was  orthogonal  to  two  others,  which 
enable  high  resolution  photography  of  the  plasma  plume  during  operation,  as  shown  in 
Figure  16. 


Figure  16.  Voltage-current  measurement  test  setup. 


With  the  necessary  infrastructure  in  place,  V-I  curves  were  then  produced  by 
setting  the  anode  flow  rate  to  a  constant  value  and  adjusting  the  discharge  voltage 
between  140  and  250  V.  At  each  voltage  the  associated  discharge  current  was  noted  and 
tabulated.  This  procedure  was  repeated  for  five  different  flow  rates.  In  specific  cases 
where  the  discharge  current  did  not  settle  to  a  single  stable  value,  six  different  values 
would  be  recorded  in  order  to  capture  the  range  of  currents  displayed. 
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Faraday  Probe 

The  Faraday  probe  was  used  in  this  analysis  to  measure  the  current  density  in  the 
plasma  plume  generated  by  the  thruster.  The  probe  was  a  straightforward  device  which 
consisted  of  a  collector  disk  that  had  a  potential  applied  to  it  and  a  resistor  that  was  used 
to  measure  the  drop  in  voltage  (F)  across  it.  The  collector  plate  was  biased  positive  in 
order  to  repel  low  energy  charge-exchange  ions  existing  in  the  plume  that  should  not  be 
included  in  the  measured  discharge  voltage.  The  Faraday  probe  used  in  this  experiment 
was  a  guarded  design  that  included  an  enclosure  surrounding  the  collector  plate  that  was 
biased  negative  in  order  to  repel  electrons  and  an  aperture  with  a  diameter  of  0.762  mm, 
in  this  case,  to  allow  ions  to  enter  and  strike  the  collector  disk,  which  is  illustrated  in 
Figure  17.  Thus,  the  following  formula  was  used  to  determine  the  current  density  (Jg) 
[43]. 

.  _  F-IOOO  /'mA\ 

^'^apeture  \CTn^) 

Where  R  is  the  value  of  the  resistor  being  used,  4.66  MOhm  in  this  case,  and  A^peture  is 
the  area  of  the  probe  guard  aperture.  The  collector  plate  was  biased  to  +30  V  using  a 
Kiethley  6517A  electrometer/high  resistance  meter  and  the  probe  body  was  biased  to  -30 
V  using  an  Agilent  6018A  power  supply.  The  voltage  drop  across  the  probe’s  resistor 
was  measured  using  an  Agilent  34970A  data  acquisition  unit  accurate  to  ±0.02%  and 
transmitted  to  a  LabView®  VI  designed  to  record  the  data  to  a  tab  delineated  file. 
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a)  Picture  of  guarded  Faraday  probe  used  in  this  test.  b)  Schematic,  from  Rotter  [42]. 

Figure  17.  Faraday  probe  depiction  and  schematic. 

Special  care  was  taken  to  achieve  accurate  alignment  of  the  thruster  with  the 
Faraday  probe.  The  method  used  was  to  first  center  the  Faraday  probe  along  the  R  stage, 
then  to  move  the  probe  all  the  way  forward  with  respect  to  the  Z  stage,  and  finally  to  set 
the  0  stage  to  0,  which  oriented  the  thruster  along  the  length  of  the  chamber.  A 
DeWalt^M  laser  level  was  then  mounted  to  a  Manfrotto'^'^  tripod  and  oriented  such  that 
the  laser  dot  was  centered  on  the  back  of  the  Faraday  probe.  The  Faraday  probe  was  then 
moved  to  the  far  corner  and  the  thruster  moved  along  the  Aerotech’s®  Y-  and  Z-axes 
until  the  laser  dot  was  centered  in  the  thruster  channel.  Finally,  the  laser  level  was 
removed  and  the  distance  between  the  Faraday  probe  aperture  and  thruster  face  precisely 
measured  with  calipers,  so  the  thruster  could  be  moved  forward  to  a  position  of  Z  =  0 
with  respect  to  the  Faraday  probe  stages.  The  first  three  steps  in  this  process  are  detailed 
below  in  Figure  18. 
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a)  Laser  level  positioned. 


b)  Laser  centered. 


c)  Thruster  alignment. 


Figure  18.  Thruster  alignment  for  Faraday  probe  data  collection. 
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The  Faraday  probe  was  attaehed  to  a  translation  system  that  enabled  the  Faraday 
probe  to  traverse  a  semieircle  around  the  thruster  about  the  plane  biseeting  the  thruster 
along  its  vertieal  axis.  Measurements  were  taken  every  1°  along  constant  radii  of  13  and 
15  cm,  as  illustrated  in  Figure  19.  The  software  automatically  calculated  the  0  value 
necessary  to  keep  the  Faraday  aperture  pointed  at  the  center  of  the  thruster  channel  exit, 
which  was  also  the  center  of  rotation.  The  software  interface  for  controlling  the  Faraday 
probe  and  collecting  the  associated  data  is  captured  in  Figure  20. 


Figure  19.  Schematic  of  Faraday  probe  experimental  setup  with  illustration  of  13  and  15 

cm  sweep  radii. 
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Figure  20.  Faraday  probe  control  and  collection  LabView  interface. 


Manipulation  of  the  Faraday  data  yielded  a  current  density  value  at  each  point  of 
measurement,  which  could  then  be  integrated  numerically  to  yield  an  estimate  for  the 
total  current  through  a  theoretical  half-sphere  shell  of  a  given  sweep  radius.  The 
assumption  was  made  that  half  of  the  hemisphere  would  be  characterized  by  the  sweep 
values  left  of  the  thruster  centerline  and  half  of  the  hemisphere  would  be  characterized  by 
the  sweep  values  right  of  the  thruster  centerline.  Thus,  the  total  beam  current  was  given 
by 

4  =  Hs^hdS  (^)  (15) 

■^0 

Additionally,  the  authors  found  it  necessary  to  characterize  the  effective 
asymmetric  bias  of  the  beam  current,  which  for  the  purpose  of  this  work  was  denoted  as 
(i^).  This  new  term  was  used  to  capture  the  portion  of  the  divergent  beam  current  that 
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would  not  cancel  based  on  a  comparison  of  the  integration  of  the  left  half-hemisphere  and 
right  half-hemisphere,  respectively,  as  shown  in  (16). 


Fa  = 
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Faraday  probe  error  was  eomputed  by  eonducting  five  conseeutive  repeat  runs  at 
identieal  test  radii,  typieally  13  and  15  em.  An  average  beam  eurrent  for  each  run,  at 
each  test  radius,  was  then  ealculated.  The  standard  deviation,  a,  was  then  eomputed  for 
the  averages  at  eaeh  radii.  A  conservative  estimate  for  the  measurement  error  assoeiated 
with  the  Faraday  probe  was  produeed  by  eomputing  the  standard  error,  as  illustrated  in 
(17). 

SE  =  t-F-  (17) 

yn 

A  t  value  of  1 .96  was  used  in  order  to  produce  a  eonfidence  factor  of  95%.  In  all  cases, 
the  value  for  n  was  5. 

Generally,  alignment  error  would  be  determined  using  quantification  of  the  full 
width  half  max  value  (FWHM)  and  a  comparison  made  between  where  this  eurrent 
density  value  was  achieved  on  one  side  of  the  plume  versus  the  other.  In  this  ease,  the 
extreme  asymmetry  present  in  plumes  made  this  teehnique  unreliable.  Extreme  care  was 
taken  in  eaeh  case  to  achieve  consistent  and  aceurate  alignment  of  the  thruster  using  the 
laser  level  teehnique  described  above.  Based  on  previous  attempts  to  use  this  method 
with  other  thrusters,  it  is  estimated  that  the  shift  was  not  more  than  two  degrees  in  any 
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case.  This  level  of  misalignment  is  well  within  the  angular  tolerance  of  the  Faraday 
probe  [44], 

Inverted  Pendulum  Stand 

AFIT’s  inverted  pendulum  stand  was  purpose  built  by  Busek  for  eleetric 
propulsion  thrust  measurement.  The  eore  of  this  instrument  was  a  linear  variable 
differential  transformer  (LVDT)  that  was  used  to  measure  the  linear  motion  of  the 
thruster  support  stand  in  response  to  an  applied  force.  It  was  discovered  that  the  initial 
resting  position  of  the  thrust  stand  placed  the  LVDT  well  outside  of  its  linear  response 
range,  so  the  LVDT  eollar  had  to  be  manually  adjusted  in  order  to  get  the  LVDT  within  a 
range  where  the  fine  adjustment  screw  eould  be  used  to  zero  the  instrument.  This  task 
required  disassembly  of  the  thrust  stand  and  repositioning  of  the  collar,  which  was 
secured  by  a  set  screw.  The  new  position  for  the  LVDT  eaused  some  concern,  beeause  it 
neeessitated  that  the  LVDT  collar  was  almost  completely  protruding  from  its  mounting 
bracket,  as  illustrated  in  Figure  21.  However,  no  apparent  issues  arose  as  a  result  of  this 
placement.  After  the  collar  was  relocated,  the  fine  adjustment  serew  was  used  to  ensure 
the  instrument  was  centered  within  its  linear  range  between  each  test. 
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Figure  21.  LVDT  collar  protruding  from  circular  mounting  bracket. 


The  output  from  the  LVDT  was  transferred  via  a  bundled  cable  and  cannon  plug 
chamber  pass  through  to  a  readout  unit  developed  by  Busek,  which  is  shown  in  Figure 
22.  The  unit  included  a  Schaevitz  MP2000  LVDT  readout  controller  that  was  used  in 
multimeter  mode  [45].  The  amplified  analog  voltage  output  from  the  MP2000  was 
recorded  by  the  ‘Thrust  Stand’  LabView  software  interface.  The  stand  also  included  a 
passive  eddy-current  damper  that  was  used  when  mounting  the  thruster  in  order  to  reduce 
large  oscillations.  The  damper  did  not  quell  the  small  random  vibrations  that  were 
transmitted  through  the  chamber,  so  it  was  not  engaged  during  calibration  runs  or  thrust 
measurement.  Another  element  on  the  thrust  stand  control  unit  was  an  inclinometer 
readout  that  provided  information  from  the  inclinometer  mounted  to  the  thrust  stand’s 
lower  platen.  The  inclinometer  was  used  to  provide  a  feedback  circuit  when  adjusting  the 
stand  level  with  a  stepper  motor,  and  the  inclination  was  noted  during  all  tests.  A 
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previous  effort  using  this  thrust  stand  revealed  that  changes  in  the  stand  level  produced  a 
highly  correlated  bias  that  could  be  corrected  if  this  data  was  noted  during  testing.  The 
LVDT  would  often  drift  slightly  during  testing,  but  an  attempt  was  made  to  minimize  this 
drift  through  the  use  of  active  water  cooling. 


Figure  22.  Thrust  stand  readout  unit. 


As  discussed  in  chapter  2,  the  repeatability  of  the  data  produced  by  the  thrust 
stand  was  highly  dependent  on  an  assumption  of  constant  temperature.  The  stand 
included  the  liberal  use  of  water  cooling  lines  on  an  external  thermal  jacket  and 
surrounding  thruster  pedestal,  as  illustrated  in  Figures  23  and  24.  The  external  water 
jacket  was  designed  to  eliminate  radiative  heat  transfer  from  the  thruster  plume  and  other 
thermal  sources  in  the  chamber.  The  lines  surrounding  the  thruster  pedestal  were  meant 
to  eliminate  any  conductive  heat  transfer  that  might  pass  from  the  thruster  body  to  the 
mounting  plate.  Cooling  water  was  supplied  from  a  NESLAB  RTE7  water  chiller  outside 
the  chamber,  which  provided  a  temperature  stability  of  ±0.01  degrees  Celsius  [46]. 
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Figure  23.  Thrust  stand  with  water  jacket  and  wiring  harness  in  place. 


Figure  24.  Thrust  stand  internal  workings. 
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The  thrust  stand  was  installed  in  the  vaeuum  chamber  on  a  large  steel  platen  in 
order  to  provide  a  stable  base  for  performing  thrust  measurements,  as  illustrated  in  Figure 
25.  The  cathode  was  mounted  on  an  external  sting  in  order  to  reduce  the  complexity  of 
thruster  integration.  After  some  experimentation  with  the  collection  of  calibration  curves 
it  was  discovered  that  adding  additional  weight  to  the  thrust  stand  platen  helped  to  reduce 
the  magnitude  of  the  oscillations  resulting  from  random  noise  transferred  through  the 
chamber  to  the  thrust  stand.  For  testing,  two  0.5  kg  weights  were  included  with  the 
thruster  on  the  mounting  pedestal. 


Figure  25.  Thrust  stand  installed  in  chamber  with  BHT-20  mounted  on  pedestal  and 
cathode  mounted  on  an  external  sting  to  reduce  complexity. 

The  thruster  was  calibrated  before  and  after  each  attempt  to  measure  thrust.  The 

process  involved  using  a  stepper  motor  that  was  controlled  via  a  switch  on  the  thruster 
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readout  unit  to  load  and  unload  weights  suspended  on  fishing  line.  These  weights  passed 
over  a  large,  low  friction,  pulley  in  order  to  translate  their  vertical  motion  to  a  horizontal 
applied  force  on  the  thrust  pedestal.  The  combined  weight  of  each  weight  plus  the 
associated  length  of  fishing  line  was  0.134  grams,  which  provided  calibration  over  a  0- 
5.2  mN  range.  The  weights  used  were  specifically  selected  for  this  test  in  order  to  reduce 
the  magnitude  of  individual  calibration  steps  when  compared  to  previous  setups.  The 
thinnest  fiexures  and  lightest  return  spring  available  for  the  thruster  stand  were  used  in 
order  to  maximize  the  resolution  of  each  calibration  step.  This  produced  reasonable 
calibration  curves,  as  evidenced  by  Figure  26;  however,  the  noise  fioor  inherent  in  this 
instrument  ultimately  proved  too  large  to  collect  any  useable  thrust  data,  as  shown  in 
Figure  27. 


Figure  26.  Thrust  stand  calibration  curve  showing  some  hysteresis,  as  illustrated  by 

Lab  View  interface. 
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Figure  27.  Thrust  stand  random  noise  eollected  over  20  minutes. 

Torsion  Balance  System 

AFIT’s  torsion  balanee  system  was  designed  and  built  by  Busek  to  measure  the 
low  thrust  levels  generated  by  micro  pulsed  plasma  thrusters  and  similar  devices.  The 
author  turned  to  this  device  after  it  became  apparent  that  the  inverted  pendulum  stand 
simply  did  not  possess  the  resolution  to  measure  the  sub-mN  forces  produced  by  the 
BHT-20.  Several  design  features  of  the  torsion  balance  collate  to  produce  an  extremely 
sensitive  and  stable  thrust  measurement  platform.  At  the  heart  of  this  instrument  is  a 
Philtec  model  D125-AB  fiber  optic  laser  displacement  sensor  (LDS),  as  shown  in  Figure 
28. 
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Figure  28.  Laser  displacement  sensor  mounted  on  micrometer  translation  stage. 


The  LDS  had  two  linear  ranges  to  select  between  for  the  desired  application  [47]. 
The  linear  range  used  in  this  study  had  a  sensitivity  of  -15.3  mV/mil.  The  other  range 
was  so  sensitive- 17  times  the  response  rate  of  the  range  selected — that  it  was  virtually 
useless  given  the  magnitude  of  the  low  frequency  oscillations  introduced  by  vibration 
through  the  chamber.  The  LDS  proved  to  be  incredibly  stable  and  was  especially  useful, 
because  it  provided  a  method  to  account  for  displacement  drift  over  time.  This  was 
essential  in  order  to  account  for  what  the  current  distance  between  the  two  calibration 
electrodes  was.  This  was  important  because  the  attractive  force,  Fcai,  generated  by  the 
electrodes  was  a  function  of  the  potential  difference  induced  between  the  electrodes,  Vcai, 
and  the  distance,  d,  between  the  electrodes,  as  shown  in  (18)  [47]. 


^cal 


(1.262XlO“^)(l+138.lXd)Fcai^ 


(18) 
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A  view  of  calibration  electrode  design  is  provided  in  Figure  29.  The  distance  between 
these  calibration  electrodes  was  adjustable  in  10  pm  increments  using  a  micrometer 
translation  stage  attached  to  the  smaller  electrode. 


Figure  29.  Torsion  balance  calibration  electrodes. 


Calibrations  runs  were  performed  immediately  before  and  after  thrust 
measurements  in  order  to  establish  the  baseline  position  and  response  function.  A 
calibration  voltage  of  800  V  was  used  in  all  cases.  It  was  generated  by  cycling  power 
(3.082  V  with  a  2  A  limit)  from  an  Agilent  E363 1 A  power  supply  that  was  routed 
through  an  EMCO  model  E121  high  voltage  amplifier  to  step  the  voltage  to  800  V.  The 
initial  distance  between  the  electrodes  was  set  to  1.5  mm  ±5  pm,  and  the  initial  EDS 
voltage  was  noted.  Changes  in  EDS  voltage  were  directly  attributable  to  changes  in 
position,  so  a  corrected  distance  was  calculated  as  part  of  each  calibration  run.  The 
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analog  voltage  from  the  LDS  was  read  using  an  Agilent  54622D,  and  the  voltages  were 
recorded  using  Agilent’s  IntuiLink  Data  Capture  software,  as  illustrated  in  Figure  30. 
Data  collection  periods  of  50  s  were  used  for  calibration  runs,  and  individual  thrust 
measurements. 


Figure  30.  Torsion  balance  data  acquisition  setup. 


Special  care  had  to  be  taken  when  integrating  the  BHT-20  onto  the  torsion 
balance,  because  of  the  instrument’s  extreme  sensitivity.  The  center  of  mass  of  the 
thruster  was  determined  by  balancing  each  axis  on  a  knife  edge.  The  COM  of  the 
thruster  was  then  located  above  the  center  of  thruster  platen  on  one  side  of  the  balance 
arm,  and  a  counterbalance  was  placed  equidistance  on  the  opposite  side  of  the  balance 
arm.  The  counterbalance  was  made  by  stacking  washers  above  a  vacuum  pass-through 
cover  until  the  combined  total  weight  was  within  2.5  grams  of  the  thruster,  as  shown  in 
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Figure  31.  A  large  single  loop  of  FTTP  propellant  line  was  used  to  transfer  xenon  to  the 
anode  in  order  to  produce  a  stable  spring  constant,  as  shown  in  Figure  32.  The  cathode 
was  mounted  externally  to  the  torsion  balance,  as  illustrated  by  Figure  33. 


Figure  31.  Thruster  counterbalance. 


Figure  32.  Back  view  of  thruster  integration  onto  torsion  balance. 
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Figure  33.  Front  view  of  thruster  integration  onto  torsion  balance. 


Error  Analysis 

The  first  element  of  error  that  was  quantified  for  the  torsion  balance  was  the 
standard  error  associated  with  the  thrust  measurements.  At  each  operating  point,  five  50 
s  collections  were  performed  at  40  Hz.  Each  of  these  data  collections  was  then  averaged 
and  the  standard  deviation,  a,  associated  with  these  averages  calculated.  The  standard 
error,  with  a  confidence  factor  of  95%,  was  then  calculated  using  the  method  discussed 
previously. 

The  calibration  bias,  Abcai,  was  determined  for  each  operating  point.  A  range  of 
two  to  seven  calibration  runs  were  performed  for  each  operating  point,  and  it  was 
desirable  to  capture  the  variation  in  resulting  calibration  slopes  values,  yi.  This  was  done 
by  first  calculating  the  average  of  the  calibration  slopes,  y,  and  then  using  (19)  to 
calculate  Abcai. 
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The  error  assoeiated  with  the  thruster  drift,  Abd,  was  also  quantified.  The 
calculation  of  Abd  consisted  of  a  comparison  of  the  voltage  baseline  before,  Vzero-,  and 
after  thrust  measurements,  Vzero+,  using  (20). 


The  error  for  the  resolution  bias,  Abres,  of  the  LDS  was  also  calculated.  The 
Agilent  data  capture  software  captured  voltage  to  1  mV.  The  resolution  bias  was 
obtained  by  dividing  that  value  in  half  and  then  multiplying  by  the  measured  change  in 
voltage  produced  by  firing  the  thruster.  In  all  cases,  this  change  in  voltage  was  less  than 
the  change  in  voltage  produced  by  the  calibration  force. 

The  total  thrust  stand  measurement  error  was  calculated  using  the  relationship 
expressed  in  (21). 


error  = 


(21) 
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IV.  Results 


This  chapter  enumerates  the  performanee  results  of  the  BHT-20  derived  from  the 
experiments  outlined  above.  The  voltage-eurrent  eharaeteristics  of  this  thruster  are 
outlined.  A  maximum  suggested  operating  power  is  provided  for  the  thruster.  Faraday 
data  is  reviewed  for  multiple  eases  in  order  to  make  attributions  about  beam  eurrent, 
beam  effieieney,  and  the  influence  of  eathode  position  on  plume  divergenee.  Thrust 
measurement  results  are  tabulated  and  reported  both  with  and  without  the  inclusion  of  the 
BHC-1500  in  order  to  differentiate  the  influenee  of  the  cathode  on  the  results,  espeeially 
with  relation  to  kp  and  effieieney. 

Voltage-Current  Characteristics 

The  anode  eurrent  and  voltage  were  measured  at  five  flow  rates  during  multiple 
thruster  runs.  In  all  of  these  cases  the  cathode  and  anode  were  both  set  to  a  flow  rate  of  1 
seem  on  startup  and  the  diseharge  voltage  was  initially  set  to  170  V.  Onee  thruster 
ignition  oeeurred  the  anode  flow  rate  was  reduced  and  the  current  allowed  to  stabilize. 
The  diseharge  voltage  was  then  varied  from  148  to  248  V,  while  the  anode  flow  rate  was 
held  at  a  constant  flow  rate.  This  proeedure  was  repeated  for  five  flow  rates.  As 
expeeted,  the  rate  at  whieh  diseharge  eurrent  inereased  with  respeet  to  an  inerease  in 
anode  potential  became  more  pronouneed  as  the  flow  rate  inereased.  It  was  also 
observed  that  at  some  operating  points  the  thruster  did  not  settle  to  a  stable  or  repeatable 
value.  This  was  the  case,  for  example,  when  the  thruster  was  operating  at  0.7  seem  and 
180  V,  as  illustrated  in  Figure  34.  It  was  believed  that  the  eause  of  the  fluetuations  was 
due  to  ionization  instability  at  those  operating  points,  but  further  testing  would  be 


55 


required  to  confirm  this  point.  Notably,  this  behavior  seemed  to  occur  over  a  limited 
discharge  voltage  range.  The  plot  reveals  that  the  BHT-20  provided  many  stable 
operating  points  spanning  the  range  of  5  -  25  Watts. 


Figure  34.  Voltage-current  characteristics  of  BHT-20. 

The  experimental  setup  used  for  the  collection  of  voltage-current  data  also 
enabled  the  simultaneous  capture  of  a  high  resolution  image  of  the  thruster  at  each 
operating  condition.  These  side-view  images  illustrate  how  the  plasma  plume  shape 
became  more  pronounced  and  extended  farther  from  the  thruster  face  as  the  discharge 
voltage  was  increased  from  148  to  248  V  at  the  constant  anode  flow  rate  of  0.45  seem. 
The  two  extreme  cases  are  shown  in  Figure  35  and  36.  The  remaining  nine  images 
associated  with  the  0.45  seem  curve  are  included  in  Appendix  B.  The  appreciable  pink 
glow  discharge  in  these  images  was  a  byproduct  of  the  comparatively  significant  1  seem 
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flow  rate  through  the  cathode  that  was  being  used  during  testing.  The  typical  blue  glow 
discharge  associated  with  ionized  xenon  is  difficult  to  separate  at  the  low  voltage  cases, 
but  it  becomes  increasingly  pronounces  as  the  discharge  voltage  was  increased.  In  all 
cases  the  discharge  looked  reminiscent  of  a  ‘ball  plume’  mode  for  a  traditional  Hall  effect 
thruster.  As  noted,  the  central  portion  of  the  plume  shape  is  difficult  to  identify  at  low 
voltages,  but  the  shape  high  density  discharge  beam  was  readily  apparent  at  higher 
discharge  voltages.  Further  insight  into  the  reason  for  the  plume  shape  seen  in  these 
images  is  included  with  the  Faraday  results  section.  An  attempt  was  initially  made  to 
extend  the  voltage  current  curves  beyond  the  range  presented  in  Figure  33,  but  it  was 
found  that  this  was  highly  detrimental  to  thruster  longevity. 


Figure  35.  BHT-20  plume  with  Vd  of  138  V  and  ma  of  0.45  seem. 
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Figure  36.  BHT-20  plume  with  Vd  of  248  V  and  ma  of  0.45  seem. 


BHT-20  Operating  Limits 

The  BHT-20  was  originally  designed  with  the  intention  to  operate  with  a 
discharge  power  of  20  W.  Early  experimentation  was  done  to  determine  the  effect  of 
operating  beyond  this  point  in  order  to  understand  if  there  would  be  any  implications  for 
thruster  longevity.  It  was  quickly  discovered  that  operating  at  25  Watts  for  hours 
produced  negligible  degradation  of  the  thruster’s  anode  post,  while  operating  at  45  Watts 
rapidly  degraded  the  post.  After  three  hours  of  operation  with  power  levels  ranging  from 
20-45  Watts,  the  BN  anode  post  was  eroded  to  the  point  where  the  mounting  screw  that 
held  it  in  place  was  mostly  exposed,  as  illustrated  by  Figure  37.  If  thruster  operation 
were  to  continue  at  this  power  level  for  an  extended  period  of  time  then  it  is  expected  that 
anode  damage  would  occur  and  thruster  performance  would  degrade.  Fortunately,  the 
anode  post  was  designed  to  be  easily  replaced  in  order  to  facilitate  such  experimentation. 
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a)  Thruster  operating  at  peak  of  45  W  b)  Resulting  damage  to  BN  anode  post. 

Figure  37.  Damage  to  thruster  after  operation  beyond  design  limit. 

Alternatively,  when  the  thruster  was  operated  over  a  range  of  5-25  Watts  for  more 
than  70  hours,  no  apparent  degradation  in  the  condition  of  the  anode  post  occurred,  as 
shown  in  Figure  38.  It  is  proposed  that  this  is  an  area  of  interest  that  should  be  monitored 
in  all  future  testing  due  to  the  necessity  of  low  power  electric  propulsion  devices  to 
compensate  for  low  thrust  levels  with  reasonable  Isp  values  and  long  operational 
lifetimes.  For  obvious  reasons,  further  testing  was  limited  to  less  than  25  W. 


a)  Thruster  operating  at  peak  of  25  W  b)  No  damage  to  BN  anode  post. 


Figure  38.  No  damage  to  thruster  after  extended  operation  at  design  limit. 
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Faraday  Data 

The  Faraday  probe  is  a  well-established  instrument  that  has  been  widely  used  to 
characterize  electric  thrusters.  The  two  most  important  traditional  values  that  can  be 
derived  from  Faraday  data  are  the  beam  current  and  divergence  angle.  The  beam  current 
can  be  used  to  calculate  a  beam  efficiency  that  relates  how  effectively  the  thruster  is 
using  the  current  flowing  from  the  cathode  to  anode  to  generate  the  ions  that  produce  the 
plasma  plume.  The  divergence  angle  has  important  implications  for  spacecraft 
integration,  especially  with  regard  to  components  that  are  especially  susceptible  to 
sputtering  erosion,  such  as  solar  panels.  The  divergence  angle  is  also  related  to  the  thrust 
factor  that  conveys  how  efficiently  momentum  transferred  to  ions  in  the  plume  will 
translate  to  propulsive  thrust. 

The  experimental  techniques  for  conducting  Faraday  probe  measurements  of  a 
Hall  effect  thruster  are  well  established,  which  is  always  ideal  when  trying  to  establish 
baseline  performance  values  for  a  new  thruster  design.  For  this  study,  Faraday  probe 
measurements  were  taken  at  different  power  levels,  thruster  orientations,  and  with  the 
cathode  in  different  positions  relative  to  the  thruster  in  order  to  identify  how  these 
changes  affected  the  performance  measures  described  above.  Further,  the  shapes  of  the 
current  density  versus  view  angle  plots  are  useful  for  making  comparisons  with  existing 
Hall  effect  thruster  designs.  However,  before  useful  comparative  Faraday  data  was 
collected,  it  was  decided  that  an  appraisal  of  the  noise  floor  should  be  performed. 

The  first  Faraday  probe  data  collected  was  a  measurement  of  the  background 
noise  level  in  the  chamber.  This  was  done  at  a  vacuum  level  of  9.7  •  10'  torr  without 
current  or  propellant  flowing  to  the  thruster  and  cathode.  The  derived  average  current 
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densities  at  radii  of  13  and  15  cm  were  -2.05  •  10'^  and  -2.02  •  10'^  respectively. 

These  averages  were  more  than  three  orders  of  magnitude  below  the  peak  current  density 
measured  for  the  thruster  when  operating  at  5  Watts,  and  the  two  distributions  were 
random  without  any  statistical  correlation  to  distance,  as  illustrated  by  Figure  39.  Given 
these  factors,  no  offset  correction  was  performed  on  the  measured  data  to  correct  for  the 
noise  floor.  This  was  a  conservative  omission  as  a  correction  to  the  data  would  have 
resulted  in  a  positive  offset.  Once  an  understanding  of  the  noise  inherent  in  the 
instrument  was  obtained,  an  examination  of  appropriate  test  radii  was  conducted. 
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Figure  39.  Background  noise  for  Faraday  probe  measurements. 

The  second  Faraday  test  was  used  to  select  and  validate  two  sweep  radii  that  were 
utilized  throughout  the  remainder  of  Faraday  measurements  in  order  to  ease  direct 
comparisons  of  current  density  distributions  at  various  operating  conditions.  The  current 
density  in  a  plasma  plume  reduces  quadratically  with  increased  distance  from  the  source, 
which  necessitated  the  selection  of  measurement  radii  that  produced  reasonable  signal  to 
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background  noise  ratios.  The  low  power  output  of  the  BHT-20  called  for  performing 
Faraday  sweeps  at  a  redueed  distanee  from  the  thruster  face  when  compared  to  previous 
researeh  performed  at  AFIT  using  the  higher  power  BHT-200  and  BHT-1500  thrusters. 
However,  previous  researeh  performed  using  this  Faraday  probe  also  suggested  that  if 
measurements  were  taken  too  elosely  to  the  source  then  the  thruster  itself  would 
influence  the  measured  data  due  to  the  generation  of  electromagnetic  noise  [42].  Initial 
Faraday  runs  were  taken  at  constant  radii  varying  between  30  and  40  cm,  but  the  resulting 
eurrent  levels  were  deemed  lower  than  was  desirable.  Additional  testing  at  radii  of  13 
and  15  em  yielded  aceeptable  measured  current  levels  for  all  of  the  operating  points  of 
interest. 

The  derived  values  generated  from  the  average  of  five  successive  Faraday  sweeps 
at  13  and  15  em  were  then  eompared  with  the  derived  values  generated  from  five 
suecessive  sweeps  at  1 1  and  17  em.  The  objective  was  to  determine  if  the  eurrent  density 
distribution,  beam  efficiency,  and  divergence  angle  remained  eonsistent  between  sweeps 
that  were  condueted  at  different  times  and  different  radii.  The  results  from  this  test 
proved  to  be  sufficiently  consistent  to  instill  confidence  in  the  validity  of  using  13  and  15 
em  radii  for  the  remainder  of  Faraday  measurements,  as  illustrated  by  Table  2  and  Figure 
40.  The  beam  effieiency,  pb,  dropped  slightly  due  to  plume  loss  factors  and  divergence 
angle,  0,  increased  slightly  due  to  ion  repulsion.  The  magnitude  of  the  current  density 
plot  deelined  appropriately  with  distance,  as  substantiated  by  the  consisteney  between  the 
integrated  beam  current  values,  Ib,  and  the  shape  of  the  distribution  remained  consistent 
as  the  measurement  radius  inereased.  For  comparison,  the  value  for  the  beam  efficiency 
was  approximately  40%  aeross  all  cases,  but  an  efficiency  of  approximately  60%  would 
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have  been  expected  for  a  highly  tuned  200  Watt  class  annular  Hall  effect  thruster  [48]. 
Further,  traditional  Hall  thrusters  tend  to  produce  highly  symmetric  discharge  plumes,  so 
it  was  decided  that  the  incongruencies  between  the  left  and  right  halves  of  Figure  39 
warranted  further  investigation. 

Table  2.  Derived  beam  parameters  for  4  W  operating  condition. 


Vacuum 

(torr) 

Radius 

(cm) 

Power 

(W) 

Vd 

(V) 

Id 

(mA) 

Ib 

(mA) 

fib 

(%) 

0 

(deg) 

Ea 

SE 

(%) 

9.51E-07 

11 

4.05 

168.7 

24 

9.77 

40.71 

69.16 

0.052 

0.38 

9.32E-07 

13 

4.05 

168.7 

24 

9.68 

40.32 

68.92 

0.056 

0.24 

9.32E-07 

15 

4.05 

168.7 

24 

9.63 

40.11 

68.95 

0.048 

0.08 

9.51E-07 

17 

4.05 

168.7 

24 

9.63 

40.10 

70.01 
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Figure  40.  Current  density  curves  for  4  W  operating  condition,  showing  consistent 
profdes  and  appropriate  quadratic  falloff 
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The  current  density  profiles  generated  from  operation  at  4  W  illustrated  two 
initially  surprising  attributes.  The  first  was  that  the  peak  current  density  measurement  did 
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not  occur  along  the  centerline  of  the  thruster,  as  was  typieal  of  annular  Hall  effect  thruster 
designs  [49].  The  “halo”  shape  was  also  not  typieal  of  curricular  Hall  effect  thruster 
designs  with  eleetromagnetic  coils,  which  also  generated  their  peak  current  density  along 
the  thruster  centerline,  but  with  a  higher  divergenee  angle  than  was  typical  of  annular 
designs  [50].  Instead,  the  distribution  was  reminiscent  of  the  results  generated  by  circular 
Hall  effect  thrusters  with  permanent  magnets,  whieh  have  current  density  profiles 
characterized  by  twin  peaks  with  a  substantially  lower  centerline  value  [51].  Extensive 
studies  of  the  latter  phenomenon  by  the  Prineeton  Plasma  Physics  Laboratory  suggest  that 
the  cause  of  the  halo  shape  resulted  from  inconsistencies  in  the  magnetic  field  topology  in 
the  vicinity  of  the  discharge  ehannel  exit  [52].  It  was  believed  that  the  same  explanation 
was  also  a  viable  cause  of  the  relative  asymmetry  between  the  left  and  right  halves  of  the 
eurrent  density  distribution. 

Additional  Faraday  measurements  were  performed  at  4.7  and  13.6  W  operating 
points  in  order  to  determine  if  the  halo  distribution  and  asymmetry  were  consistently 
present.  The  thrust  asymmetry  faetor.  Fa,  for  these  cases  is  provided  in  Table  3  and  4. 

The  results  suggest  that  thrust  asymmetry  and  the  general  shape  observed  at  4  W  were 
consistent  issues,  but  the  degree  of  asymmetry  present  in  the  13.6  W  case  was  less  than 
half  that  of  the  4  or  4.7  W  ease.  Unfortunately,  the  beam  efficiency  and  divergence  angle 
were  both  less  desirable  in  the  13.6  W  case.  A  graph  of  the  current  density  distribution 
for  4.7  and  13.6  W  operating  points  is  provided  in  Figure  41  and  42.  The  shapes  of  the  4 
and  4.7  W  distributions  were  remarkably  similar,  whereas  the  shape  of  the  right  peak  for 
the  13.6  W  case  was  substantially  more  uniform-without  any  appreciable  step  on  the 
right  half  of  the  plume. 
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Table  3.  Derived  beam  parameters  for  4.7  W  operating  condition. 


Vacuum  Radius  Power  Vd  Id  Ib  ilb  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa  (%) 

l.OOE-06  13  4.72  168.7  28  11.79  42.09  69.24  0.058  2.48 

l.OOE-06  15  4.72  168.7  28  11.70  41.78  69.43  0.053  1.55 


Table  4.  Derived  beam  parameters  for  13.6  W  operating  condition. 


Vacuum  Radius  Power  Vd  Id  Ib  Ob  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa  (%) 

1.25E-06  13  13.57  178.6  76  27.17  35.75  73.92  0.032  0.77 

1.25E-06  15  13.57  178.6  76  27.14  35.71  74.20  0.029  0.44 


Angle  Off  Centerline  (degrees) 

Figure  41.  Current  density  curves  for  4.7  W  operating  condition. 
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Angle  Olf  Centerline  (degrees) 

Figure  42.  Current  density  eurves  for  13.6  W  operating  eondition. 


Conelusions  about  eurrent  density  asymmetry  from  the  first  three  eases  had  one 
significant  limitation.  The  measurements  were  all  taken  along  the  same  bisecting  plane 
and  the  asymmetry  phenomenon  was  inherently  a  three  dimensional  attribute.  Thus,  in 
order  to  gain  a  better  understanding  of  the  asymmetry,  the  4.7  and  13.6  W  cases  were 
repeated  after  the  thruster  had  been  rotated  90  degrees  clockwise,  with  respect  to  the  front 
face.  The  thrust  asymmetry  factor  for  these  cases  is  provided  in  Table  5  and  6.  The 
results  for  both  of  these  cases  were  consistently  lower  thrust  asymmetry  factors  than  for 
the  values  reported  above  at  the  same  respective  power  levels.  This  suggests  that  many 
rotations  at  finer  degree  increments  would  be  required  in  order  to  accurately  define  the 
true  thrust  asymmetry.  Similarly,  the  beam  efficiency  improved  for  the  4.7  W  case  and 
degraded  for  the  13.6  W  case,  with  respect  to  the  previous  measurements,  and  the 
divergence  angle  broadened  for  the  4.7  W  case  and  narrowed  slightly  for  the  13.6  W 
case.  These  observations  provide  further  motivation  for  future  plume  studies  on  this 
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thruster  to  include  numerous  thruster  rotations.  A  graph  of  the  current  density 
distribution  for  these  4.7  and  13.6  W  operating  points  is  provided  in  Figure  43  and  44. 
The  general  shape  of  these  distributions  was  consistent  with  the  previous  test,  but  fine 
details  near  the  peaks  illustrate  subtle  changes  that  resulted  from  the  90  degree  thruster 
rotation.  The  distribution  for  the  rotated  4.7  W  case  illustrates  a  broadening  of  the  left 
peak  and  an  increase  in  the  right  plateau  relative  to  the  right  peak.  The  distribution  for 
the  rotated  13.6  W  case  illustrates  a  five  degree  shift  towards  the  centerline  for  the  left 
peak.  These  changes  serve  to  reinforce  the  validity  of  the  changes  in  the  reported  Fa 
values.  It  was  suggested  that  the  cathode  position  may  have  some  influence  on  the  shape 
of  the  current  density  curves,  so  a  study  was  devised  to  look  at  the  impact  of  varying 
cathode  position  in  small  increments. 

Table  5.  Derived  beam  parameters  at  4.7  W  operating  condition  with  thruster  rotated  90 
degrees  with  respect  to  configuration  for  Table  3. 


Vacuum  Radius  Power  Vd  Id  Ib  hb  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa  (%) 

9.73E-07  13  4.72  168.7  28  1 1. 14  39.78  67.61  0.023  O.IO 

9.73E-07  15  4.72  168.7  28  1 1. 13  39.77  67.75  0.020  0.09 


Table  6.  Derived  beam  parameters  at  13.6  W  operating  condition  with  thruster  rotated  90 
degrees  with  respect  to  configuration  for  Table  4. 


Vacuum 

Radius 

Power 

Vd 

Id 

Ib 

hb 

0 

SE 

(torr) 

(cm) 

(W) 

(V) 

(mA) 

(mA) 

(%) 

(deg) 

Ea 

(%) 

1.20E-06 

13 

13.75 

178.6 

77 

27.75 

36.03 

75.05 

0.006 

0.19 

1.20E-06 

15 

13.75 

178.6 

77 

27.77 

36.06 

75.12 

0.004 

0.14 
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Angle  Off  Centerline  (degrees) 

Figure  43.  Current  density  curves  for  4.7  W  operating  condition  with  thruster  rotated 
90  degrees  with  respect  to  configuration  for  Figure  40. 


Angle  Off  Centerline  (degrees) 


Figure  44.  Current  density  curves  for  13.6  W  operating  condition  with  thruster  rotated 
90  degrees  with  respect  to  configuration  for  Figure  41. 

Variation  of  Cathode  Position 

The  BHT-20  was  initially  operated  with  the  hollow  cathode  limited  to  a  single 
position;  the  cathode  exhaust  port  was  1.59  cm  above  the  centerline  of  the  thruster 
discharge  channel  and  1.59  cm  from  the  front  face  of  the  thruster.  This  placement  would 
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pose  particular  challenges  for  integration  onto  a  very  small  spacecraft,  so  it  was  desirable 
to  look  at  the  impact  of  moving  the  cathode  relative  to  the  thruster.  The  cathode  was 
moved  farther  away  from  the  thruster  centerline  in  1  mm  increments  over  a  range  of  1.59 
to  2.09  cm.  The  cathode  was  also  moved  relative  to  the  thruster  face  in  1  mm  increments 
over  a  range  of  1 .39  to  1 .89  cm.  The  two  axes  the  cathode  was  moved  along  are 
illustrated  in  Figure  45.  In  all  cases  the  anode  mass  flow  rate  was  held  constant  at  0.48 
seem. 


Figure  45.  Axes  used  for  movement  of  the  cathode  relative  to  the  thruster  face  and 

centerline. 

It  was  found  that  the  shape  of  the  current  density  plot  and  divergence  angle  were 
the  only  two  derived  performance  measures  that  significantly  changed  over  the  1 1  test 
cases  conducted  at  varying  cathode  positions.  Results  from  the  two  extreme  points  are 
included  below  in  Table  7  and  8  and  Figure  46  and  47  to  illustrate  the  point.  It  is  notable 
that  increasing  the  distance  from  the  centerline  consistently  increased  the  divergence 
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angle  and  decreasing  the  distance  relative  to  the  thruster  face  consistently  increased  the 
divergence  angle,  as  illustrated  by  Figure  48.  Thus,  the  divergence  angle  reported  in 
Table  7  shows  a  nearly  six  degree  reduction  in  divergence  angle  when  compared  to  the 
values  in  Table  8,  It  was  believed  that  changes  in  the  cathode  coupling  voltage  were 
responsible  for  the  observed  relationship  between  cathode  position  and  divergence  angle 
based  on  previous  results  from  Sommerville  and  King  on  Hall  effect  thrusters  [53].  The 
left  peak  of  the  current  density  plot  for  Figure  44  was  shifted  nearly  10  degrees  closer  to 
the  thruster  centerline  relative  Figure  45.  The  plateau  in  Figure  44  became  the  right  peak, 
whereas  the  right  peak  for  Figure  45  was  consistent  with  previous  results  at  4.7  W. 
Quantitative  results  from  the  additional  test  cases  are  included  in  Appendix  C.  The  next 
step  in  the  characterization  of  the  thruster  was  to  expand  the  Faraday  data  to  capture  the 
extreme  points  of  the  BHT-20  operating  range. 

Table  7.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  1.59  cm 
above  thruster  centerline  and  1 .29  cm  from  thruster  face. 

Vacuum  Radius  Power  Vd  Id  Ib  ^h  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa  (%) 

9.39E-07  13  4.72  168.7  28  11.20  39.99  64.99  0.021  1.02 

9.39E-07  15  4.72  168.7  28  11.20  39.99  64.56  0.016  0.93 


Table  8.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  2.09  cm 
above  thruster  centerline  and  1.59  cm  from  thruster  face. 

Vacuum  Radius  Power  Vd  Id  Ib  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Ea  (%) 

9.39E-07  13  4.72  168.7  28  11.06  39.52  70.07  0.021  0.82 

9.39E-07  15  4.72  168.7  28  1 1. 14  39.77  70.94  0.022  0.97 
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Current  Density  (mA/ern^)  31  Current  Density  (mA/ern^) 


Angle  Otf  Centerline  (degrees) 

46.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  1.59  cm 
above  thruster  centerline  and  1 .29  cm  from  thruster  face. 


Angle  Off  Centerline  (degrees) 

Figure  47.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  2.09  cm 
above  thruster  centerline  and  1.59  cm  from  thruster  face. 
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Figure  48.  Plot  of  divergence  angle  as  a  function  of  distance  from  the  X  and  Y  axis, 
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As  evidenced  by  the  previous  test  cases,  the  author  had  a  strong  interest  in  the 
ability  of  the  BFlT-20  to  operate  at  power  levels  below  the  rated  20  W  for  the  thruster. 
This  concentration  stemmed  from  the  suggestion  that  the  thruster  was  a  possible 
candidate  for  very  small  satellites  with  highly  constrained  power  budgets.  Nonetheless,  it 
was  also  prudent  collect  Faraday  data  over  the  entire  thruster  operating  range  in  order  to 
understand  how  the  derived  performance  parameters  varied  with  discharge  power. 

On  the  low  end  of  the  scale,  it  was  possible  to  stably  and  repeatedly  operate  the 
thruster  down  to  a  power  level  of  2.5  W.  The  derived  results  for  that  operating  condition 
are  provided  in  Table  9  and  Figure  49.  The  values  presented  are  consistent  with  those 
reported  for  operation  at  4  and  4.7  W,  but  it  was  noteworthy  that  the  beam  efficiency  was 
the  highest  of  any  operating  condition  tested  over  the  entirety  of  this  study.  Similarly,  the 
shape  of  the  current  density  plot  was  similar  to  the  4  and  4.7  W  cases,  but  a  noticeable 
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narrowing  of  the  beam  was  visible.  On  the  opposite  end  of  the  speetrum,  a  20  W  ease 
was  identified  for  eomparison  purposes. 

Table  9.  Derived  beam  parameters  for  2.5  W  operating  eondition. 


Vacuum 

Radius 

Power 

Vd 

Id 

Ib 

hb 

0 

SE 

(torr) 

(cm) 

(W) 

(V) 

(mA) 

(mA) 

(%) 

(deg) 

Ea 

(%) 

l.OOE-06 

13 

2.53 

148.7 

17 

7.26 

42.71 

69.70 

0.035 

0.53 

l.OOE-06 

15 

2.53 

148.7 

17 

7.31 

43.03 

69.87 

0.032 

0.62 

Angle  Off  Centerline  (degrees) 

Figure  49.  Current  density  eurves  for  2.5  W  operating  condition. 

The  channel  size  and  SmCo  magnet  were  originally  selected  based  on  a  20  W 
operating  condition,  so  it  was  anticipated  that  the  derived  performance  parameters  would 
peak  around  that  operating  condition.  Flowever,  it  was  found  that  the  only  performance 
parameter  that  was  optimal  at  20  W  was  the  thrust  asymmetry  term,  which  reached  a 
minimum  value  of  0.3%.  That  was  roughly  one  tenth  the  degree  of  asymmetry  found  for 
the  2.5  W  operating  case.  Unfortunately,  the  beam  efficiency  and  divergence  angle  were 
only  marginally  better  than  the  13.6  W  case,  but  worse  than  all  others.  The  results  are 


73 


provided  in  Table  10  and  Figure  50.  The  current  density  distribution  was  similar  to  the 
shape  observed  at  13.7  W,  but  the  distribution  appeared  to  broaden  and  a  double  peak 
manifested  itself  on  both  sides  of  the  20  W  case.  It  was  observed  during  testing  of  the 
extreme  cases  that  it  was  possible  to  get  the  thruster  to  operate  in  an  additional — more 
diffuse — mode,  so  the  performance  associated  with  the  new  mode  was  measured. 

Table  10.  Derived  beam  parameters  for  20  W  operating  condition. 

Vacuum  Radius  Power  Vd  Id  Ib  fib  0  SE 

(torr)  (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa  (%) 

I.OOE-06  13  20.02  238.3  84  31.00  36.90  74.35  0.006  0.26 

I.OOE-06  15  20.02  238.3  84  30.90  36.79  74.64  0.003  1.22 


Angle  Off  Centerline  (degrees) 


Figure  50.  Current  density  curves  for  20  W  operating  condition. 

Comparison  of  Plume  Modes 

During  testing  the  author  discovered  it  was  possible,  with  careful  control  of  the 
anode  flow  rate,  to  settle  the  thruster  to  the  same  discharge  power  level  with  different 
mass  flow  rates  through  the  anode.  This  was  generally  achieved  after  thruster  startup  by 
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setting  the  discharge  voltage  to  a  desired  level  and  then  crashing  the  mass  flow  rate  from 
1  seem  to  a  flow  rate  slightly  below  that  which  was  generally  necessary  to  achieve  a 
desired  discharge  current  when  incrementally  decreasing  the  flow  rate.  In  all  of  the  4.7 
W  cases  that  were  detailed  above,  the  mass  flow  rate  was  approximately  0.45  seem.  The 
derived  values  detailed  in  Table  1 1  and  Figure  51  were  for  a  4.7  W  operating  condition 
where  the  anode  flow  rate  was  only  4  seem. 

Table  11.  Derived  beam  parameters  for  4.7  W  operating  condition  with  reduced  anode 

mass  flow  rate. 


Vacuum 

Radius 

Power 

Vd 

Id 

Ib 

fib 

0 

SE 

(torr) 

(cm) 

(W) 

(V) 

(mA) 

(mA) 

(%) 

(deg) 

Fa 

(%) 

9.42E-07 

13 

4.72 

168.7 

28 

6.91 

24.66 

80.90 

0.053 

0.18 

9.42E-07 

15 

4.72 

168.7 

28 

6.78 

24.23 

80.71 

0.057 

0.26 

Angle  Off  Centerline  (degrees) 


Figure  5 1 .  Current  density  curves  for  4.7  W  operating  condition  with  reduced  anode 

mass  flow  rate. 

The  ability  to  achieve  the  same  discharge  power  with  a  lower  mass  flow  rate  is 
always  desirable,  but  in  this  case,  the  cost  of  reducing  the  flow  rate  was  significant 


75 


degradation  of  the  other  performance  measures.  The  shape  of  the  current  density  plot 
revealed  significantly  reduced  peak  currents  when  compared  to  higher  flow  rate  cases  and 
the  centerline  value  was  reduced  to  the  level  of  the  extreme  wings.  The  result  of  the 
reduction  in  measured  current  density  was  a  significantly  lower  estimate  for  the  total 
beam  current.  In-turn,  this  reduced  the  beam  current  efficiency  by  approximately  40% 
when  compared  to  the  higher  flow  cases.  The  other  significant  attribute  was  a 
broadening  of  the  plasma  plume,  which  was  illustrated  by  the  largest  divergence  angle 
measured  for  any  test  case.  This  test  concluded  the  Faraday  analysis  performed  on  the 
thruster. 

Thrust  Measurement 

Measurement  of  the  thrust  produced  by  the  BHT-20  was  especially  challenging  due 
to  the  relatively  low  magnitude  of  the  force  generated  by  the  thruster  when  operating  over 
a  range  of  2.8  to  21 .9  W.  An  attempt  was  made  to  measure  thrust  with  the  inverted 
pendulum  stand  detailed  above,  but  the  noise  inherent  in  the  instrument  simply  dwarfed 
the  signal.  This  necessitated  reconfiguring  the  chamber  to  accept  the  torsion  balance, 
which  proved  to  be  the  ideal  instrument  for  this  application.  Four  different  successful 
test  cases  were  conducted  and  their  results  summarized. 

The  first  experiment  performed  with  the  torsion  balance  involved  conducting  14 
calibration  runs  with  the  cathode  on  and  the  thruster  off.  The  cathode  was  left  on  to 
ensure  that  its  operation  did  not  influence  the  linearity  assumption  for  the  LDS  response 
to  a  change  in  force.  The  14  cases  consisted  of  two  runs  at  seven  different  calibration 
force  conditions.  The  seven  cases  were  produced  by  setting  the  electric  potential  between 
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the  calibration  electrodes  values  of  200  to  800  V,  in  100  V  increments.  When  plotted,  the 
result  was  a  curve  of  LDS  voltage  change  as  a  function  of  applied  calibration  force,  as 
depicted  in  Figure  52.  The  curve  that  resulted  had  a  least  squares  regression  value  of 
0.9996  and  spanned  beyond  the  range  of  measured  thrust  values  for  the  BHT-20,  which 
provided  positive  validation  of  resolution  and  linearity  of  torsion  balance.  Therefore,  it 
was  decided  that  thrust  testing  could  commence  with  the  torsion  balance. 
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Figure  52.  LDS  voltage  change  as  a  function  of  applied  calibration  force. 

The  results  provided  in  this  section  disregarded  the  mass  flow  rate  associated  with 
the  BHC-1500  cathode,  because  a  much  different  cathode  would  be  paired  with  the  BHT- 
20  for  operational  applications.  Therefore,  the  value  for  Igp  included  below  is 
systematically  higher  than  was  actually  observed  during  this  research  effort,  but  it  serves 
to  provide  a  reasonable  estimate  for  the  hp  that  would  result  if  the  thruster  were  paired 
with  a  thermionic  cathode.  Results  from  all  previous  attempts  to  scale  down  Hall 
thrusters  to  low  power  levels  suggested  that  efficiency  would  decrease  nonlinearly,  so  it 
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was  not  surprising  that  thrust  data  from  the  BHT-20  validated  that  general  rule,  as  shown 
in  Table  12. 

Table  12.  BHT-20  measured  thrust  and  derived  values  using  the  torsion  balance. 


Anode  Potential 
(V) 

Anode  Current 
(A) 

Power 

(W) 

Mass  Flow 
(Pg/s) 

Thrust 

(pN) 

Error 

(%) 

ha 

(%) 

Isp 

(s) 

148.73 

0.019 

2.83 

44.24 

44 

3.22 

0.76 

100.72 

168.62 

0.028 

4.72 

44.24 

75 

2.23 

1.36 

173.58 

178.56 

0.120 

21.43 

60.95 

170 

1.00 

1.23 

299.42 

238.50 

0.092 

21.94 

49.15 

181 

3.77 

1.52 

375.27 

The  anode  efficiencies  reported  for  the  BHT-20  are  lower  than  desirable  and 
suggest  further  optimization  is  necessary  prior  to  operational  use.  Nonetheless,  the  hp 
values  do  offer  some  promise  that  with  optimization  the  thruster  could  be  a  viable 
contender  for  some  missions.  It  was  noteworthy  that  the  highest  hp  was  achieved  at  the 
highest  discharge  voltage,  which  is  consistent  with  electrostatic  thruster  theory.  Circular 
Hall  effect  thrusters  also  tend  to  perform  better  at  relatively  high  discharge  voltages, 
often  above  300  V,  so  this  suggests  further  testing  should  be  performed  at  higher 
voltages. 

Impact  of  Cathode  on  Performance 

The  results  provided  in  Table  13  illustrate  why  it  would  not  be  reasonable  to  pair 
the  BHT-20  with  a  cathode  that  was  originally  designed  to  be  paired  with  a  1500  W 
annular  Hall  effect  thruster.  In  several  cases  the  required  xenon  flow  rate  for  the  cathode 
was  more  than  twice  the  flow  being  used  by  the  anode  to  generate  ions,  which  is 
evidenced  by  the  cathode  efficiency  values,  pc.  Accounting  for  the  total  mass  flow  means 
that  the  actual  Isp  values  were  31.26,  52.08,  1 14.59,  and  125.09  s  for  the  2.83,  4.72, 
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21.43,  and  21.94  W  operating  cases,  respectively.  To  make  matters  worse,  the  electrical 
efficiency,  po,  was  also  extremely  poor  due  to  the  amount  of  power  required  by  the 
cathode.  These  two  factors  were  compounded  when  the  total  efficiency,  px,  was 
calculated. 

Table  13.  BHT-20  measured  thrust  and  derived  values  using  the  torsion  balance. 


Cathode  Anode 


Keeper 

Potential 

(V) 

Keeper 

Current 

(A) 

Cathode 

Power 

(W) 

Discharge 

Power 

(W) 

Mass 

Flow 

(Pg/s) 

Mass 

Flow 

(Pg/s) 

he 

(%) 

ho 

(%) 

hr 

(%) 

21.5 

0.6 

12.90 

2.83 

98.30 

44.24 

31.04 

17.97 

0.04 

21.3 

0.6 

12.78 

4.72 

98.30 

44.24 

31.04 

26.98 

0.11 

20.5 

0.6 

12.29 

21.43 

98.30 

60.95 

38.27 

63.55 

0.27 

20.2 

0.6 

12.13 

21.94 

98.30 

49.15 

33.33 

64.40 

0.33 

Findings 

The  BHT-20  represented  a  step  forward  in  the  ability  to  operate  electrostatic 
thrusters  at  very  low  discharge  powers.  The  thruster  was  found  to  have  many  stable 
operating  points  over  the  range  of  2.5  to  25  W,  and  no  degradation  in  the  thruster  was 
observed  while  operating  at  those  power  levels.  However,  operation  at  45  W  caused 
rapid  erosion  of  the  anode  post,  so  the  remainder  of  tests  focused  on  operating  points  in 
the  2.5  to  20  W  region.  Plume  sweeps  conducted  with  a  Faraday  probe  revealed  that  the 
thruster  shares  many  attributes  with  small  permanent  magnet  circular  Hall  thruster 
designs.  The  plume  shape  exhibited  a  characteristic  halo  distribution  with  a  relatively 
broad  divergence  angle  when  compared  with  efficient  annular  Hall  thruster  designs.  It 
was  also  note  that  at  certain  operating  points  the  distribution  of  the  ions  was  quite 
asymmetric,  which  is  not  a  desirable  attribute,  but  it  could  likely  be  curtailed  with 
magnetic  shunting,  as  described  in  chapter  5.  It  was  observed  that  variation  of  the 
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cathode  position  over  a  narrow  range  eaused  ehanges  in  the  current  density  distribution 
profile  and  divergence  angle,  but  did  not  signifieantly  impact  any  other  performance 
measures.  Thrust  measurements  revealed  that  the  BHT-20  had  a  lower  fip  and  total 
efficieney  than  the  preponderance  of  Hall  thruster  designs  previously  employed; 
however,  to  a  degree  this  was  expected  based  on  empirical  results  with  sealing  Hall 
thrusters  to  low  power.  Of  course,  the  observed  values  were  signifieantly  influenced  by 
the  use  of  an  unnecessarily  large  and  power  hungry  cathode  for  this  applieation. 
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V.  Conclusions  and  Recommendations 


This  chapter  presents  preliminary  findings  based  on  this  body  of  work. 
Conelusions  regarding  the  BHT-20’s  performanee  are  put  forth  and  reeommendations  for 
further  work  with  this  thruster  are  provided.  Possible  avenues  for  further  improvement  of 
the  thruster  and  test  instruments  are  suggested.  The  BHT-20  was  both  a  prototype  and  a 
novel  departure  from  previous  Hall  effect  thruster  design,  so  a  broad  range  of  questions 
remain  to  be  answered  by  future  efforts. 

Nonetheless,  all  of  the  initial  objeetives  of  this  research  effort  were  met.  The 
BHT-20  was  operated  at  diseharge  powers  ranging  from  2.5  to  20  W.  Operating 
proeedures  and  integration  ehallenges  were  identified,  viable  solutions  implemented,  and 
doeumentation  generated  for  future  efforts.  Voltage-eurrent  curves  were  produeed  at  five 
flow  rates  and  stable  operating  points  identified  for  investigation  with  the  techniques 
previously  deseribed.  Faraday  measurements  were  performed  and  the  assoeiated 
derivative  values  reported.  The  ion  current  density  profiles  revealed  the  thruster  behaved 
like  a  permanent  magnet  eireular  Hall  thruster  in  terms  of  plot  shape  and  ealeulated 
divergenee  angle.  Thrust  was  measured  at  four  operating  points,  and  this  was  used  to 
generate  values  for  fip  and  effioieney.  The  values  produeed  by  this  analysis  were  lower 
than  assoeiated  with  nearly  all  Hall  effect  thrusters  designed  to  operate  at  20  to  100  W, 
but  the  measured  values  do  suggest  that  further  investigation  is  warranted  for  CubeSat 
applieations  where  most  eurrent  eleetrie  thrusters  are  simply  too  power  hungry  to  be 
viable. 
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Voltage  Current  Curve  Implications 

The  voltage-current  curves  generated  as  part  of  this  research  effort  provide  a 
useful  tool  for  further  investigations,  but  it  is  recommended  that  additional  data  be 
gathered  to  resolve  some  of  the  peculiarities  presented  in  this  study.  The  development  of 
complete  VI  curves  was  complicated  by  the  fact  that  anode  post  was  found  to  rapidly 
degrade  when  operated  near  40  W.  No  current  limited  mode  was  identified  for  the 
thruster  at  any  flow  rates,  because  the  author  did  not  desire  to  push  the  thruster  above  its 
20  W  design  point  and  kept  too  relatively  conservative  discharge  voltages.  In  retrospect, 
the  similarities  exhibited  between  the  operation  of  the  BHT-20  and  a  circular  Hall 
thruster  suggest  that  the  thruster  may  have  performed  better  at  higher  discharge  voltages. 
For  this  reason  it  is  suggested  that  at  low  flow  rates,  which  do  not  result  in  the  thruster 
exceeding  25  W,  the  VI  curves  be  expanded  by  including  discharge  voltages  up  to  400  V. 

The  second  part  of  the  VI  curves  that  deserve  more  study  are  the  unstable 
operating  points.  Data  acquisition  software  should  be  used  to  capture  the  anode  current 
fluctuations  at  those  operating  points.  The  author  suggests  initially  starting  with  a  100 
Hz  sampling  rate,  which  would  enable  the  experimenter  to  study  behavior  occurring  at 
below  50  Hz  according  to  Nyquist  criterion.  It  is  proposed  the  instability  may  be  a  result 
of  high  frequency  breathing  modes,  but  further  investigation  will  be  required  to  test  this 
hypothesis. 

BHT-20  Thruster 

This  thruster  represents  the  first  Hall  effect  variant  to  successfully  operate  in  the 
2.5-5  W  range.  The  fact  the  thruster  will  operate  stably  and  without  apparent  degradation 


82 


over  a  range  that  spans  2.5  to  25  W  is  no  less  signifieant.  The  effieieneies  and  hp  values 
achieved  over  this  range  are  lower  than  would  warrant  immediate  integration  onto  any 
operational  system,  but  with  further  refinement  it  is  possible  that  this  thruster  design  will 
evolve  into  a  viable  propulsion  source  for  very  small  power  limited  satellites.  Unless 
substantial  improvements  to  all  of  the  performance  parameters  are  made  then  it  is  not 
suggested  that  this  design  be  considered  as  a  viable  alternative  in  the  10+  W  category  due 
to  the  wide  availability  of  competitive  alternatives.  Thus,  one  logical  next-step  would  be 
to  refine  the  understanding  of  this  technology  through  the  formulation  and  execution  of  a 
design  of  experiments  aimed  at  better  characterizing  the  current  performance  envelope  of 
this  design. 

The  performance  characterization  in  this  effort  explored  the  results  for  a  wide 
range  of  operating  conditions,  but  it  was  not  a  parametric  study  intended  to  identify 
optimal  operating  points  in  terms  of  fip,  thrust,  or  efficiency.  It  is  likely  that  with  further 
investigation  operating  points  could  be  identified  that  exceed  the  best  performance  results 
found  in  this  study.  Nonetheless,  it  is  suggested  that  such  an  effort  be  limited  to  low 
power  points  in  order  to  reduce  the  trade-space  and  focus  the  study  where  this  thruster 
design  holds  the  possibility  to  mature  into  a  viable  operational  technology.  It  is 
suggested  that  before  such  an  effort  is  undertaken  a  better  suited  cathode  be  paired  with 
the  BHT-20. 

The  impact  of  the  BHC-I500  cathode  on  the  total  efficiency  and  fip  values  found 
in  this  study  reveal  that  the  only  way  this  thruster  would  be  of  any  operational  value  is  if 
it  were  paired  with  a  thermionic  cathode.  HeatWave  Labs,  Inc.  is  currently  the  author’s 
leading  candidate  to  supply  such  a  device,  but  other  manufacturers,  such  as  Kymball 
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Physics  and  Heat  Wave  Labs,  Ine.,  also  make  single  erystal  emitters  that  are  large  enough 
to  generate  an  eleetron  eurrent  of  greater  than  30  mA.  It  is  known  that  pairing  a  single 
erystal  LaBe  would  eause  nearly  all  of  the  performanee  measurements  in  this  study  to 
ehange,  but  it  is  not  currently  possible  to  aecurately  predict  exactly  how  these  measures 
would  ehange  due  to  a  laek  of  prior  work  in  this  area.  It  is  known  that  the  elimination  of 
the  1  seem  of  xenon  that  was  supplied  to  the  eathode  would  yield  significant 
improvements  in  the  kp,  cahode  effieieney,  and  total  efficieney  terms.  It  is  also  known 
that  the  power  required  to  operate  an  appropriately  sized  single  crystal  emitter  will  be 
lower,  so  nearly  all  of  the  performanee  numbers  are  expeeted  to  improve.  This,  eoupled 
with  the  BHT-20’s  ability  to  operate  at  discharge  power  levels  below  any  previous  Hall 
thruster  design  provides  a  eompelling  argument  to  pursue  further  experimentation. 

One  of  the  most  exciting  discoveries  that  oeeurred  while  operating  the  BHT-20 
was  its  ability  to  funetion  stably  at  diseharge  powers  down  to  2.5  W.  This  suggests  it 
would  be  possible,  when  paired  with  a  thermionie  cathode,  to  serve  as  a  propulsion 
system  on  a  3U  CubeSat.  The  eleetrieal  effieieney  is  eoneeming,  beeause  it  suggests  a 
signifieant  percentage  of  the  applied  power  is  going  to  loss  meehanisms,  some  of  which 
are  produeing  heating  that  may  be  difficult  to  manage  on  a  CubeSat.  If  possible,  design 
variations  should  be  pursued  that  aim  to  improve  the  eleetrieal  efficieney  of  this  thruster, 
but  an  inerease  in  the  kp  would  likely  be  seen  as  an  even  more  important  area  of  foeus. 
Small  satellites  must  be  good  stewards  of  their  available  power,  but  at  least  the  power 
eomes  from  an  external  souree.  The  kp  relates  how  effeetively  the  propellant  that  is 
brought  with  the  satellite  is  being  used  to  produee  a  momentum  exehange.  The  latter  is  an 
area  where  future  evolutions  of  this  design  eould  potentially  eompete  effeetively  against 
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the  current  competition  in  the  5  W  thruster  class.  A  continuous  operating  power  of  5W  is 
deemed  reasonable  for  integration  onto  current  CubeSats,  because  AFIT  is  currently 
preparing  to  launch  a  CubeSat  with  a  20  W  average  orbital  power  that  is  generated  via  a 
deployable  seven-panel  solar  array. 

The  BHT-20  is  rife  with  possible  areas  for  experimentation.  Increasing  the 
magnetic  field  strength,  altering  the  frit  design,  and  resizing  the  thruster  channel  around 
operation  at  5  W-instead  of  20  W-all  hold  promise.  The  wide  range  of  suggestions  is 
possible,  in  part,  due  to  the  relatively  simple  and  inexpensive  design  of  the  BHT-20.  For 
comparison,  this  thruster  could  be  modified  to  become  a  circular  Hall  thruster  with  the 
integration  of  a  magnet  in  place  of  the  BN  anode  post  at  the  center  of  the  discharge 
channel.  Of  course,  this  suggestion  would  only  make  sense  if  the  additional  magnet 
could  be  maintained  at  temperatures  well  below  its  Curie  temperature.  Alternatively,  this 
magnet  could  be  integrated  into  the  center  of  the  anode,  below  the  anode  post,  where  it 
might  be  better  protected  from  the  thermal  environment,  but  would  still  provide  an 
appreciable  radial  component  to  the  magnetic  field  lines  within  the  thruster  channel.  This 
would  enable  the  thruster  to  generate  a  magnetic  trap  via  the  electron  gyration  radius  that 
is  expected  to  increase  the  ionization  efficiency  of  the  thruster.  One  immediate 
suggestion  to  improve  the  current  design  is  to  shunt  the  magnetic  field  lines  outside  of  the 
channel  exit,  because  this  has  been  shown  to  eliminate  the  halo  plume  associated  with 
permanent  magnet  circular  Hall  thruster  designs  [51].  In-turn,  one  could  reasonably 
expect  an  improvement  in  the  thrust  factor  and  thrust  asymmetry  factor  associated  with 
the  current  design.  The  latter  could  be  verified  with  additional  Faraday  sweeps. 
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Faraday  Probe 

The  Faraday  probe  provided  useful  and  repeatable  measurements  that  were 
fundamental  in  establishing  a  preliminary  theory  for  how  the  BHT-20  functioned.  One  of 
the  important  insights  provided  by  Faraday  probe  measurements  was  the  asymmetry 
typically  associated  with  most  operating  points.  However,  the  Fa  term  was  derived  based 
on  measurements  from  a  single  plane  passing  through  the  centerline  of  the  thruster,  and 
measurements  performed  along  two  different  planes  displayed  different  Fa  values,  as  was 
expected.  It  is  suggested  a  rotisserie  be  constructed  which  could  rotate  the  thruster  about 
its  centerline  without  the  need  to  open  the  chamber  between  measurements.  This  would 
enable  the  rapid  generation  of  a  three-dimensional  estimate  for  Fa.  Ultimately,  this  could 
be  used  to  predict  the  magnitude  and  direction  of  the  actual  thrust  asymmetry  inherent  in 
this  design.  Similarly,  the  other  instruments  used  in  this  characterization  would  also 
benefit  from  some  minor  modifications. 

Inverted  Pendulum  Thrust  Stand 

The  inverted  pendulum  thrust  stand  proved  to  be  a  reliable  instrument,  but  the 
author  suggests  that  its  use  be  limited  to  the  measurement  of  thrust  levels  above  2  mN  in 
order  to  achieve  a  reasonable  signal  to  noise  ratio.  One  suggested  improvement  for  the 
current  thrust  stand  configuration  is  the  integration  of  a  higher  flow  water  pump  in  order 
to  keep  the  heat  rejection  rate  constant;  the  current  setup  was  unable  to  keep  a  constant 
flow  rate  over  the  timespan  necessary  to  complete  a  single  calibration  run.  Another 
suggestion  for  future  researchers  is  to  use  the  HD  webcam  in  the  chamber  to  view  the 
calibration  weight  placement  in  order  to  ensure  repeatable  calibration  curves.  The  author 
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had  improved  success  with  the  introduction  of  the  latter  technique.  However,  the  author 
does  not  feel  any  improvements  to  the  inverted  pendulum  stand  would  make  it  a  viable 
alternative  to  the  torsion  balance  for  the  direct  measurement  of  the  thrust  output  of  the 
BHT-20. 

Torsion  Balance 

The  torsion  balance  proved  to  be  the  ideal  tool  for  measuring  the  relatively  low 
thrust  levels  associated  with  the  BHT-20.  The  results  were  repeatable  and  the  associated 
error  was  very  low  considering  the  associated  signal.  Nonetheless,  a  couple  of  minor 
changes  could  dramatically  improve  the  user  experience.  Surfaces  within  large  vacuum 
chambers  can  change  orientation  slightly  after  the  chamber  achieves  its  operating 
pressure  and  thousands  of  pounds  of  force  are  exerted  against  its  external  structure.  This 
behavior  was  observed  with  the  inverted  pendulum  stand  via  its  inclinometer.  Such  a 
device  could  easily  be  integrated  onto  the  torsion  stand.  Of  course,  in  order  for  such  a 
tool  to  be  useful  the  two  level  adjustment  micrometers  on  the  torsion  balance  would  need 
to  have  externally  controlled  stepper  motors  attached  to  them.  Along  the  same  line,  it 
would  also  be  very  useful  to  attach  similar  stepper  motors  to  the  LDS  and  calibration 
electrode  micrometer  translation  stages.  Future  testing  of  the  BHT-20  at  AFIT  will 
undoubtedly  involve  the  inverted  pendulum  stand,  but  some  additional  instruments 
should  also  be  considered. 

Other  Diagnostic  Techniques 

This  study  was  limited  in  scope  to  the  use  of  a  Faraday  probe  and  torsion  balance, 
but  other  tools  would  be  required  in  order  to  validate  the  life  expectancy  of  the  BHT-20 
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and  to  develop  a  better  understanding  of  how  it  operates.  The  author  believes  one 
pressing  study  that  remains  is  the  thermal  eharacterization  of  this  thruster.  It  is 
speculated  that  thermal  heating  is  not  an  issue  for  this  thruster  due  to  the  fact  that  no 
discernable  change  in  magnetic  field  strength  occurred  after  more  than  70  hours  of 
testing,  but  this  fact  would  need  to  be  verified  prior  to  a  long  endurance  space  mission. 
Further,  it  would  be  desirable  to  identify  the  steady  state  temperature  of  the  anode  in 
order  to  determine  if  this  design  could  be  morphed  into  a  circular  or  end  Hall  thruster 
design  without  special  cooling  requirements.  Laser  induced  fluorescence  is  another 
technique  that  could  be  used  to  validate  the  operating  assumptions  surrounding  how  this 
thruster  generates  plasma  and  where  the  ion  acceleration  occurs.  The  theory  underlying 
both  of  these  techniques  is  included  in  Appendix  D. 


Appendix  A,  BHT-20  Operating  Procedures 

The  following  procedures  were  used  for  the  operation  of  the  BHT-20  plasma  thruster  and 
the  BHC-1500  cathode  being  used  with  it.  Please  see  operating  manuals  on  SPASS 
laboratory  computers  for  instructions  on  how  to  operate  various  test  instruments 
identified  in  Capt  de  La  Harpe’s  thesis. 

1.  Ensure  vacuum  chamber  pressure  is  below  5  •  10'^  torr. 

2.  Purge  xenon  propellant  lines; 

a.  Open  xenon  gas  bottle  valve. 

b.  Set  propellant  regulator  valve  to  20  psi. 

c.  Turn  on  MKS  four-channel  readout. 

d.  Set  flow  controllers  for  anode  and  cathode  to  1  seem. 

e.  Open  propellant  Swagelok  valves  (if  applicable), 
f  Allow  propellant  to  flow  for  30  minutes. 

g.  Secure  flow  controllers. 

3.  Condition  Cathode 

a.  If  the  cathode  has  been  exposed  to  ambient  atmospheric  conditions  for  a 
considerable  amount  of  time  (>4  days),  then  set  flow  controller  for 
cathode  to  0.5  seem  during  conditioning  process.  Otherwise  no  flow. 

b.  Turn  on  the  PPU  power  supply  and  cooling  fan. 

c.  Start  BPU-600  Lab  View  program. 

d.  Set  heater  current  to  2  A  for  90  minutes. 

e.  Increase  heater  current  to  4  A  for  90  minutes, 
f  Increase  heater  current  to  6  A  for  30  minutes. 
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g.  Turn  off  heater  current. 

4.  Prepare  Propellant  Lines 

a.  Open  mass  flow  controller  to  cathode  and  set  to  a  flow  rate  of  1  seem. 

5.  Light  Cathode 

a.  Set  cathode  heater  current  to  6.5  A  until  breakdown  voltage  reaches  5.9  V 
(~10  minutes). 

b.  Set  cathode  keeper  current  to  0.6  A. 

i.  Cathode  ignition  is  indicated  by  a  drop  in  keeper  voltage. 

ii.  If  keeper  voltage  remains  above  100  V,  cathode  failed  to  break 
down.  Repeat  from  step  5. 

c.  Secure  cathode  heater:  set  cathode  heater  to  0  A. 

6.  Light  Thruster 

a.  Set  anode  mass  flow  controller  to  1  seem. 

b.  Set  discharge  voltage  to  170  V. 

i.  Blue  glow  discharges  at  the  anode  face  indicates  the  thruster  is  lit. 

ii.  If  thruster  fails  to  light  then  confirm  that  cathode  is  lit  and 
reattempt  step  6. 

7.  Set  to  desired  discharge  power. 

a.  Use  voltage  current  curves  from  Capt  de  La  Harpe’s  Master’s  thesis  to 
identify  desired  operating  point  based  on  discharge  voltage  and  discharge 
current. 
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b.  Reduce  anode  flow  rate  to  the  flow  rate  associated  with  the  desired 
operating  point.  This  should  be  done  within  20  seconds  of  lighting 
thruster  in  order  to  ensure  thruster  does  not  dwell  at  power  level  above  25 
W. 

c.  Set  discharge  voltage  to  desired  value 

i.  Allow  current  to  decrease  from  starting  value  and  stabilize  before 
adjusting  discharge  voltage. 

ii.  Move  discharge  voltage  to  desired  value  in  10  V  incriments. 

d.  Adjust  anode  flow  rate  until  desired  discharge  current  is  achieved. 

8.  Allow  thruster  to  stabilize  and  conduct  testing.  Ensure  chamber  pressure  remains 
below  -'1x10'^  torn 

9.  Securing  Thruster 

a.  Secure  power  to  the  anode  and  keeper. 

b.  Secure  xenon  flow  to  the  anode  and  cathode. 

c.  Close  propellant  Swagelok  valves  (if  applicable). 

d.  Close  xenon  bottle  valve. 

e.  Close  Host  Simulator  software  and  secure  power  to  the  PPU. 

f  Turn  off  PPU  cooling  fans  after  allowing  PPU  to  cool  for  10  minutes, 
g.  Turn  off  all  power  supplies. 

10.  If  restoring  the  vacuum  chamber  atmospheric  pressure: 

a.  Ensure  the  cathode  has  had  at  least  60  minutes  to  cool  after  operation. 

b.  Push  red  button  on  vacuum  control  panel. 
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c.  Fill  chamber  with  dry  nitrogen  onee  eryo-pumps  have  eeased  operation. 
Ensure  that  ehamber  door  elamps  are  not  seeured  tightly  to  avoid  ereation 
of  pressure  vessel. 
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Appendix  B,  BHT-20  Plume  Pictures  at  Increasing  Discharge  Voltages, 

A  series  of  photos  are  presented  that  illustrate  how  the  diseharge  plume  changes 
as  the  discharge  voltage  is  increased  and  the  anode  flow  rate  is  held  constant  at  0.45  seem 
xenon.  The  photos  presented  here  are  associated  with  voltages  ranging  from  158  to  238 
V,  in  10  V  increments.  A  discussion  of  the  trend  across  the  entire  voltage  range  is 
included  in  chapter  4,  and  additional  insight  into  the  phenomena  visible  in  these  images  is 
provided  by  the  analysis  of  Faraday  data  included  in  chapter  4.  At  a  discharge  voltage  of 
158  V  the  cathode  discharge  significantly  obscures  the  typical  blue  illumination 
associated  with  ionized  xenon.  Close  observation  reveals  a  ball-like  plume  just 
protruding  from  the  discharge  channel  where  the  ion  density  is  sufficient  to  provide  a 
nearly  white  illumination,  as  illustrated  in  Figure  53. 


Figure  53.  BHT-20  plume  with  Vd  of  158  V  and  ma  of  0.45  seem. 
Increasing  the  discharge  voltage  to  168V  does  not  substantially  change  the  visible 
behavior  illustrated  above,  but  the  central  portion  of  the  plume  does  become  slightly 
more  prominent,  as  shown  in  Figure  54. 
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Figure  54.  BHT-20  plume  with  Vd  of  168  V  and  ma  of  0.45  seem. 


When  the  discharge  voltage  is  increased  to  178  V,  the  plume  becomes  incrementally 
more  defined,  as  illustrated  by  Figure  55. 


Figure  55.  BHT-20  plume  with  Vd  of  178  V  and  ma  of  0.45  seem. 


Again,  increasing  the  discharge  voltage  to  188  V  causes  the  high  density  discharge  beam 
to  become  more  distinctly  differentiated  from  the  cathode  discharge,  as  depicted  in 
Figure  56. 
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Figure  56.  BHT-20  plume  with  Vd  of  188  V  and  ma  of  0.45  seem. 


At  a  discharge  voltage  of  198  V,  the  ball-like  nature  of  the  plasma  plume  is  readily 
apparent  and  relatively  easy  to  differentiate  from  the  cathode  discharge,  as  illustrated  by 
Figure  57.  At  this  increment,  the  impact  of  the  electric  field  lines  on  the  cathode  plume 
also  becomes  more  pronounced. 


Figure  57.  BHT-20  plume  with  Vd  of  198  V  and  ma  of  0.45  seem. 
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The  general  trend  established  above  eontinues  to  hold  with  an  inerease  in  the  diseharge 
voltage  to  208V,  as  shown  in  Figure  58. 


Figure  58.  BHT-20  plume  with  Vd  of  208  V  and  ma  of  0.45  seem. 

The  established  trend  eontinues  to  hold  with  an  inerease  in  the  diseharge  voltage  to  218V, 
as  shown  in  Figure  59. 


Figure  59.  BHT-20  plume  with  Vd  of  218  V  and  ma  of  0.45  seem. 
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The  established  trend  eontinues  to  hold  with  an  inerease  in  the  discharge  voltage  to  228V, 
as  shown  in  Figure  60. 


Figure  60.  BHT-20  plume  with  Vd  of  228  V  and  ma  of  0.45  seem. 


The  established  trend  continues  to  hold  with  an  increase  in  the  discharge  voltage  to  238V, 
as  shown  in  Figure  61 .  The  luminosity  produced  by  the  central  plume  has  increased 
substantially  over  the  nine  cases  and  the  central  plume  is  now  easily  distinguishable. 


Figure  61.  BHT-20  plume  with  Vd  of  238  V  and  ma  of  0.45  seem. 
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Appendix  C,  Faraday  Probe  Data 


The  data  included  in  this  section  illustrates  the  impact  of  moving  the  cathode 
relative  to  the  thruster.  The  cathode’s  position  was  varied  relative  to  the  thruster 
centerline  and  thruster  face  in  1  mm  increments.  Table  14  identifies  the  specific  test 
cases  included  in  this  section.  In  all  cases  the  anode  mass  flow  rate  was  held  constant.  It 
was  found  that  the  shape  of  the  current  density  plot  and  divergence  angle  were  the  only 
two  derived  performance  measures  that  significantly  changed.  It  was  found  that 
increasing  the  distance  from  the  centerline  consistently  increased  the  divergence  angle 
and  decreasing  the  distance  from  the  thruster  face  consistently  increased  the  divergence 
angle. 


Table  14.  Test  cases  included  in  Appendix  C. 


Position 

Figure 

Number 

Table 

Number 

Distance  from  Cathode  to 
Thruster  Centerline  (cm) 

Distance  from  Cathode 
Thruster  Face  (cm) 

A 

62 

15 

1.59 

1.59 

B 

63 

16 

1.59 

1.39 

C 

64 

17 

1.59 

1.49 

D 

65 

18 

1.59 

1.69 

E 

66 

19 

1.59 

1.79 

F 

67 

20 

1.69 

1.59 

G 

68 

21 

1.79 

1.59 

H 

69 

22 

1.89 

1.59 

I 

70 

23 

1.99 

1.59 
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Position  A  represents  the  baseline  run,  beeause  it  was  performed  with  the  eathode 
in  the  position  speeified  by  the  manufaeturer  of  the  BHT-20.  The  plume  shape  revealed 
in  Figure  62  aligns  well  with  the  expeetations  established  after  operating  the  thruster  at 
numerous  operating  conditions.  The  beam  efficiency,  divergence  angle,  and  thrust 
asymmetry  factor  reported  in  Table  15  also  align  well  with  the  behavior  reported  in 
chapter  4. 


Angle  Off  Centerline  (degrees) 

Figure  62.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  A. 


Table  15.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  A. 

Vacuum  Radius  Power  Va  Id  Ib  ^b  0 

(torr) (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.40E-07  13  4.72  168.7  28  11.19  39.98  67.40  0.017 

9.40E-07  15  4.72  168.7  28  11.06  39.49  66.49  0.013 


99 


The  distance  from  the  cathode  to  the  thruster  face  was  reduced  by  2  mm  for 
position  B,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  63  illustrates  an 
appreciable  decrease  in  the  peak  current  density  measurement  and  an  appreciable 
broadening  of  the  curve.  These  visual  cues  are  substantiated  by  the  results  presented  in 
Table  16,  which  reports  that  the  divergence  angle  increased  while  the  other  figures  of 
merit  remained  relatively  constant. 


Angle  Off  Centerline  (degrees) 

Figure  63.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  B. 


Table  16.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  B. 

Vacuum  Radius  Power  Vd  Id  Ib  ^b  0 

(torr) (cm)  (W)  (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.09E-07  13  4.72  168.7  28  11.21  40.05  69.84  0.017 

9.09E-07  15  4.72  168.7  28  11.13  39.76  69.53  0.012 
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The  distance  from  the  cathode  to  the  thruster  face  was  reduced  by  1  mm  for 
position  C,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  64  illustrates  a 
slight  decrease  in  the  peak  current  density  measurement  and  a  slight  broadening  of  the 
curve.  These  visual  cues  are  substantiated  by  the  results  presented  in  Table  17,  which 
reports  that  the  divergence  angle  increased  while  the  other  figures  of  merit  remained 
relatively  constant,  with  respect  to  the  baseline  case.  The  divergence  angle  for  position  C 
was  between  those  measured  at  position  A  and  position  B,  so  this  serves  to  form  the 
initial  basis  of  a  trend. 


Angle  Off  Centerline  (degrees) 

Figure  64.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  C 


Table  17.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  C. 

Vacuum  Radius  Power  Vd  Id  Ib  fib  0 

(torr) (cm)  (W) (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.56E-07  13  4.72  168.7  28  11.17  39.88  69.00  0.017 

9.56E-07  15  4.72  168.7  28  11.16  39.87  68.83  0.014 
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The  distance  from  the  cathode  to  the  thruster  face  was  increased  by  1  mm  for 
position  D,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  65  illustrates  a 
slight  increase  in  the  peak  current  density  measurement  and  a  slight  narrowing  of  the 
curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated  by  the  results 
presented  in  Table  18,  which  reports  that  the  divergence  angle  decreased  while  the  other 
figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline. 


Angle  Off  Centerline  (degrees) 

Figure  65.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  D. 


Table  18.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  D. 


Vacuum  Radius  Power 

(torr) (cm) (W) 

9.73E-07  13  4.72 

9.73E-07  15  4.72 


Vd  Id  Ib  fib  0 

(V)  (mA)  (mA)  (%)  (deg)  Ea 

168.7  28  11.08  39.56  67.06  0.020 

168.7  28  11.07  39.53  66.83  0.015 
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The  distance  from  the  cathode  to  the  thruster  face  was  increased  by  2  mm  for 
position  E,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  66  illustrates  an 
appreciable  increase  in  the  peak  current  density  measurement  and  an  appreciable 
narrowing  of  the  curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated 
by  the  results  presented  in  Table  19,  which  reports  that  the  divergence  angle  decreased 
while  the  other  figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline. 
The  results  for  divergence  angle  as  a  function  of  distance  from  the  thruster  face  showed  a 
consistent  improvement,  over  the  range  tested,  as  the  distance  increased. 


Angle  Off  Centerline  (degrees) 

Figure  66.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  E. 


Table  19.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  E. 

Vacuum  Radius  Power  Vd  Id  Ib  fib  0 

(torr)  (cm)  (W) (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.39E-07  13  4.72  168.7  28  11.16  39.87  66.38  0.019 

9.39E-07  15  4.72  168.7  28  11.17  39.90  65.94  0.012 
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The  distance  from  the  cathode  to  the  thruster  centerline  was  increased  by  1  mm 
for  position  F,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  67  illustrates  a 
slight  decrease  in  the  peak  current  density  measurement  and  a  slight  narrowing  of  the 
curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated  by  the  results 
presented  in  Table  20,  which  reports  that  the  divergence  angle  increased  while  the  other 
figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline. 


Angle  Off  Centerline  (degrees) 

Figure  67.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  F. 


Table  20.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  F. 

Vacuum  Radius  Power  Va  Id  Ib  ^b  0 

(torr) (cm)  (W) (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.33E-07  13  4.72  168.7  28  11.06  39.51  68.88  0.021 

9.33E-07  15  4.72  168.7  28  11.08  39.56  68.51  0.017 
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The  distance  from  the  cathode  to  the  thruster  centerline  was  increased  by  2  mm 
for  position  G,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  68  illustrates 
an  appreciable  decrease  in  the  peak  current  density  measurement  and  an  appreciable 
narrowing  of  the  curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated 
by  the  results  presented  in  Table  21,  which  reports  that  the  divergence  angle  increased 
while  the  other  figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline 


Angle  Off  Centerline  (degrees) 

Figure  68.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  G. 


Table  21.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  G. 

Vacuum  Radius  Power  Vd  Id  Ib  fib  0 

(torr) (cm)  (W) (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.42E-07  13  4.72  168.7  28  11.07  39.55  69.22  0.022 

9.42E-07  15  4.72  168.7  28  11.06  39.51  68.99  0.019 
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The  distance  from  the  cathode  to  the  thruster  centerline  was  increased  by  3  mm 
for  position  H,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  69  illustrates 
an  appreciable  decrease  in  the  peak  current  density  measurement  and  an  appreciable 
narrowing  of  the  curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated 
by  the  results  presented  in  Table  22,  which  reports  that  the  divergence  angle  increased 
while  the  other  figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline. 


Angle  Off  Centerline  (degrees) 

Figure  69.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  H. 


Table  22.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  H. 


Vacuum  Radius  Power 

(torr) (cm) (W) 

9.70E-07  13  4.72 

9.70E-07  15  4.72 


Vd  Id  Ib  fib  0 

(V)  (mA)  (mA)  (%)  (deg)  Ea 

168.7  28  11.11  39.66  69.46  0.022 

168.7  28  11.07  39.55  69.04  0.017 
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The  distance  from  the  cathode  to  the  thruster  centerline  was  increased  by  4  mm  for 
position  I,  relative  to  position  A.  The  plume  shape  revealed  in  Figure  70  illustrates  an 
appreciable  decrease  in  the  peak  current  density  measurement  and  an  appreciable 
narrowing  of  the  curve,  relative  to  the  baseline  case.  These  visual  cues  are  substantiated 
by  the  results  presented  in  Table  23,  which  reports  that  the  divergence  angle  increased 
while  the  other  figures  of  merit  remained  relatively  constant,  with  respect  to  the  baseline. 
The  results  for  divergence  angle  as  a  function  of  distance  from  the  thruster  centerline 
showed  a  consistent  improvement,  over  the  range  tested,  as  the  distance  decreased. 


Angle  Off  Centerline  (degrees) 

Figure  70.  Current  density  curves  for  4.7  W  operating  condition  with  cathode  in 

position  I. 


Table  23.  Derived  beam  parameters  for  4.7  W  operating  condition  with  cathode  in 

position  I. 


Vacuum  Radius  Power  Vd  Id  Ib  ilb  0 

(torr) (cm)  (W) (V)  (mA)  (mA)  (%)  (deg)  Fa 

9.44E-07  13  4.72  168.7  28  11.07  39.55  70.18  0.025 

9.44E-07  15  4.72  168.7  28  11.03  39.40  69.40  0.019 
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Appendix  D,  Background  and  Theory  of  Further  Diagnostic  Techniques 


Thermal  Camera  Background, 

The  next  test  suggested  for  the  BHT-20  is  thermal  characterization  using 
thermocouples  and  a  FLIR  camera.  This  test  is  necessary  to  ensure  that  the  steady  state 
operating  temperature  of  the  BHT-20  does  not  approach  the  Currie  temperature  of  the 
SmCo  permanent  magnet  being  used  in  the  design,  because  that  would  cause 
demagnetization  that  would  degrade  the  thruster’s  performance.  The  infrared  camera  that 
is  suggested  is  the  ThermaCAM  SC640,  which  has  an  operational  range  of  -40°C  to 
1500°C  with  an  accuracy  of  ±2%  of  reading  [54].  This  instrument  is  illustrated  in  Figure 
71. 


Figure  71.  ThermaCAM  SC640  Infrared  Camera  [54]. 

Thermal  Camera  Theory, 

Matter  above  absolute  zero  emits  electromagnetic  radiation  with  wavelengths 
between  0.75  micrometers  and  100  micrometers  (infrared)  that  can  be  correlated  to  a 
surface  temperature  if  the  properties  of  the  matter  are  known  [4].  The  ideal  absorber  and 
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emitter  of  radiation  is  a  blaekbody.  The  radiation  emitted  by  a  black  body,  is 
idealized  and  can  be  represented  as  a  function  of  wavelength  and  temperature. 


EudT) 
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T  is  the  absolute  temperature,  h  is  Plank’s  constant,  c  is  the  speed  of  light,  and  k  is 
Boltzmann’s  constant.  Integrating  this  equation  over  all  wavelengths  produces  the 
Stefan-Boltzmann  Law  for  total  blaekbody  power. 
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cr  is  the  Stephan-Bolzmann  constant. 

KirchhofP  s  Law  relates  the  emissivity  of  a  surface  to  its  absorptivity,  where  a  is 
the  ratio  of  spectral  radiant  power  absorbed  to  incident  power  and  £  is  the  emissivity. 
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These  two  equations  can  be  combined  to  yield  a  relationship  between  the  actual  and 
measured  temperature  based  on  the  optical  characteristics  of  the  test  specimens  [55]. 
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Team  refers  to  the  blackbody  temperature  reported  by  the  thermal  eamera,  Tabj  is  the 
surface  temperature  of  the  specimen’s  surface,  Tygfi  is  the  effective  temperature  or 
radiation  incident  on  the  surface,  and  Typ  is  the  temperature  of  the  view  port.  The 
emissivity  and  transmissivity,  t,  are  a  function  of  wavelength,  temperature,  and  material. 
Emissivity  is  also  function  of  the  finish  of  the  material,  so  published  data  for  this  value 
will  need  to  be  substantiated  using  an  independent  temperature  measurement  device,  such 
a  thermocouple. 

Temperature  measurements  using  the  infrared  camera  will  require  special  care. 

The  quality  of  the  data  collected  will  be  a  product  of  alignment,  view  angle,  and  focus. 
Further,  the  vacuum  chamber’s  viewport  will  need  to  be  regularly  polished  in  order  to 
avoid  increased  thermal  noise  from  deposits  that  will  form  on  the  window  after  thruster 
operation.  Thermocouples  will  be  used  on  external  surfaces  in  order  to  have  a  source  to 
substantiate  the  temperatures  captured  by  the  infrared  camera. 

Laser  Absorption, 

Initial  visual  observation  of  the  micro-Hall  thruster  revealed  that  the  greatest 
plasma  luminescence  was  heavily  concentrated  toward  to  center  of  ionization  channel. 
This  observation  should  be  verified  using  laser  absorption.  Laser  absorption  works  based 
on  the  fundamental  principle  of  spectroscopy  that  each  ionization  state  for  a  given  atom 
or  molecule  possesses  unique  spectral  lines  which  represent  the  wavelength  of  light  that 
the  atom  or  molecule  will  absorb.  The  pattern  of  spectral  lines  observed  for  a  specific 
element  is  related  to  the  statistical  distribution  of  the  electrons  for  that  element  and  the 
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associated  quanta  of  energy,  AE,  that  will  cause  one  of  the  electrons  to  transition  from 
one  state,  Ei,  to  another,  E2,  which  is  described  by  (26)  [56]; 

AE  =  E2-  E2  (26) 


The  speeific  wavelengths  of  light  that  will  be  absorbed  are  related  to  the,  AE, 
from  (26)  through  the  theory  of  wave-particle  duality,  whieh  deseribes  the  energy  a 
photon  carries  as  a  funetion  of  its  energy  level  (E).  Thus  the  wavelength  that  will  cause  a 
speeific  electron  transition  can  be  found  using  (27)  where  h  is  Plank’s  constant,/is  the 
photon  frequency,  c  is  the  speed  of  light  and  2  is  the  wavelength  of  the  photon  [56]. 

E  =  hf=’j  (27) 


The  applieation  of  absorption  to  eharaeterizing  this  eleetric  thruster  stems  from 
the  faet  that  it  is  possible  to  couple  the  properties  above  with  some  additional 
relationships  in  order  to  determine  the  number  density  of  a  speeific  atom  within  a  given 
length.  When  a  laser  is  emitting  photons  tightly  distributed  about  the  wavelength 
assoeiated  a  strong  spectral  line  is  passed  through  a  volume  of  gas  that  absorbs  at  that 
speetral  line  the  result  is  a  decrease  in  power,  which  can  be  modeled  according  to  Beer’s 
Eaw  for  absorption,  which  is  given  as  (28)  [56].  E  is  the  beam  intensity  resulting  from 
the  energy  lost  to  absorption,  /  is  the  length  of  the  path  through  the  volume  of  absorbing 
gas,  Ni  the  number  density  of  the  atoms  and  Oabs  is  the  absorption  cross-sectional  area. 


/  =  I  ^abs^^i) 
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A21,  the  Einstein  spontaneous  emission  rate  eoeffieient,  and  gi  and  g2,  whieh  are 
degeneraey  terms  related  to  the  transitioning  level  of  the  atom  or  moleeule  in  question. 
The  ehallenge  is  that  one  needs  a  referenee  for  these  values;  however,  the  benefit  is  that 
one  ean  solve  (29)  for  number  density,  produeing  (30),  in  order  to  experimentally 
determine  the  eoncentration  of  a  speeifie  atom  or  moleeule  in  a  given  volume  by 
eomparing  the  ratio  of  the  laser  reference  intensity  with  the  reduced  intensity  measured 
after  the  laser  has  passed  through  the  gas  of  interest. 


^absl 


(30) 


The  absorption  lines  for  the  first  three  xenon  ionization  states  are  included  in 
Appendix  A.  This  study  will  focus  on  the  identification  of  the  first  ionization  state.  The 
absorbance  will  measured  by  comparing  the  transmitted  laser  intensity  to  the  etalon  using 
a  Matsui’s  methodology  [57]. 

In  order  to  take  laser  absorption  measurements  the  baseline  thruster  design  will 
need  to  be  revised.  Absorption  measurements  will  require  a  focuses  beam  of  laser  light 
to  pass  uninterrupted  through  the  xenon  plasma  within  the  ionization  channel.  The  back 
plate  of  the  thruster  is  currently  opaque  to  infrared  radiation,  so  the  design  and  material 
selection  for  the  back  plate  will  need  to  be  revised.  The  current  proposed  redesign  would 
introduce  the  neutral  xenon  flow  from  the  sides  of  the  thruster  and  include  cutouts  in  the 
back  plate  that  would  be  sealed  with  a  transparent  diamond  sheet  in  order  to  provide 
windows  for  laser  transmission.  The  fabrication  of  this  design  will  require  additional 
investigation. 
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